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EVALUATION 


The  objective  of  the  8 tody  wan  to  develop  a  method  of  determining  the 
detection  performance  of  ground-baaed  radars  against  airborne  targetH  under 
a  wide  variety  of  environment  and  system  conditions. 

This  effort  involved  the  development  of  a  novel  computationally  efficient 
method  for  simulating  ground-based  coherent  radar  system  performance 
subjected  to  target  and  clutter  signals  which  are  partially  correlated. 

Tin*  short  com  lit}*,  ol  conventional  simulation  techn  l<|iies  used  to  compute 
detection  performance  Is  that  an  extremely  large  number  of  statistical 
replications  are  required  to  establish  false  alarm  performance.  This  work 
advanced  the  art  by  adapting  the  technique  oF  importance  sampling  which 
significantly  reduces  the*  required  computer  time,  even  when  simulating 

nonlinear  systems.  In  most  cases  it  is  possible  to  simulate  false  alarm 
-8  3  4 

probabilities  as  low  as  10  with  only  about  10  to  10  replications.  A 
report,  entitled  "Importance  Sampling  Applied  to  Radar  Ka I se- Alarms, " 
summarizing  results  and  application  examples  has  been  accepted  for 
publication  in  the  IEEE  Transactions  of  the  AES. 
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1.  INTRODUCTION  AND  SUMMARY 

The  objective  of  the  study  described  in  this  report  was  to  develop  a  method 
of  determining  detection  performance  for  a  ground-based  radar  against  an  air¬ 
borne  target  under  a  wide  variety  of  environment  and  system  conditions.  This 
objective  was  met  by  constructing  a  digital  simulation  of  the  radar  signals  as 
they  are  transmitted,  reflected  from  the  scatterers  representing  targets  and 
clutter,  and  received  by  and  processed  In  the  radar  receiver.  Realistic  ground, 
rain,  and  chaff  clutter  environments  arc  modeled,  as  well  as  the  processing  fea¬ 
tures  of  modern  radars  such  as  moving  target  indication  (MT1),  both  coherent  and 
noncoherent  integration,  constant  false  alarm  rate  (CFAR)  processing,  and  non¬ 
linear  operations.  The  detection  performance  is  determined  by  Monte  Carlo 
sampling  techniques,  in  which  the  target  and  clutter  scattering  models  are 
described  statistically,  as  well  as  the  location  of  the  target  with  respect  to 
the  center  of  the  antenna  beam,  range  gate,  and  Doppler  filter. 

The  shortcoming  of  conventional  sampling  techniques  applied  to  detection 
performance  in  radar  is  the  extremely  large  number  of  statistical  replications 
required  to  establish  false  alarm  performance.  In  order  to  overcome  this  prob¬ 
lem  a  considerable  effort  was  expended  in  developing  importance  sampling  tech¬ 
niques  that  could  be  applied  to  the  wide  variety  of  signals  in  radar,  including 
the  non-Gaussian  and  non-Rayleigh  signals  characteristic  of  clutter.  The  results 
of  this  effort  are  described  in  Section  2.  In  most  cases  it  is  possible  to  simu¬ 
late  false  alarm  rates  as  low  as  10  ®  with  only  about  10^  to  10^  replications  of 
the  experiment.  Techniques  for  handling  mixed  statistics  as  well  as  CFAR  are  also 
desc  rihed . 

It  is  well  known  that  antenna  motion  during  the  coherent  processing  <.ime  of 
the  radar  causes  the  clutter  to  be  amplitude  modulated;  the  effect  is  a  broadening 
of  the  clutter  spectrum  that  would  be  observed  if  there  were  no  antenna  motion. 

It  is  straightforward  to  compute  the  resultant  spectrum  if  the  clutter  is  spa¬ 
tially  homogeneous,  especially  If  the  antenna  pattern  and  original  clutter  spec¬ 
trum  are  Gaussian  shaped.  In  Section  3  we  extend  the  analysis  to  nonhomogen eo us 
clutter  and  arbitrary  beamshapes.  While  the  signals  are  nonstationary,  the  re¬ 
sulting  algorithms  are  in  a  form  that  is  amenable  to  efficient  digital  simulation. 

In  Section  4  a  general  procedure  is  derived  for  generating  correlated  random 
signal  sequences  that  are  characteristic  of  clutter  in  a  ground-based  radar.  The 
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procedure  Is  extremely  efficient  as  it  is  based  on  fast  Fourier  transform  (FFT) 
combined  with  interpolation,  it  can  be  applied  to  any  number  of  signal  samples 
and  the  shape  of  the  power  spectral  density  is  arbitrary. 

The  simulation  program  that  was  developed  to  determine  detection  performance 
is  described  in  Section  5.  and  a  Fortran  listing  is  given  in  Appendix  A. 
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2.  IMPORTANCE  SAMPLING 


Importance  sampling  is  a  technique  that  can  be  applied  to  the  simulation 
ot  low-probability  events  without  Incurring  the  computation  costs  usually  as¬ 
sociated  with  such  simulations.  With  importance  sampling  one  can  modify  the 
probability  distribution  of  the  underlying  random  process  in  order  to  make 
the  low-probability  events  (false  alarms)  occur  more  frequently.  The  desired 
probabilities  at  the  output  of  the  process  are  then  found  by  weighting  each 
event  by  a  factor  that  is  a  function  of  only  the  state  of  the  Input;  this  fac¬ 
tor  is  independent  of  the  process  itself  [J“6]. 

The  basic  principle  of  importance  sampling  is  straightforward  as  described 
above.  However,  it  is  not  so  well  known  just  how  the  technique  can  be  made  to 
work  in  a  particular  application.  For  example,  what  if  there  are  multiple  ran¬ 
dom  inputs  to  some  processor  where  the  inputs  might  belong  to  different  statis¬ 
tical  processes?  Or  if  the  processor  is  nonlinear?  Or  if  there  is  not  a  unique 
relationship  between  the  input  and  output  of  the  processor?  These  and  other  is¬ 
sues  will  be  addressed  in  this  report .  We  will  concentrate  on  applications  to 
the  simulation  of  signals  in  radar  and  communication  systems  in  order  to  limit 
the  scope  of  the  study.  We  begin  with  a  tutorial  discussion  of  both  conventional 
and  importance  sampling. 

7 . \  CONVENTIONAL  SAMP! AHC  THEORY 

Let  us  designate  the  output  of  a  statistical  process  as  y.  We  wish  to 
estimate  the*  probability  density  function  of  this  process,  p(y),  or  its  cumu¬ 
lative  distribution  function,  P(y) ,  with  a  finite  set  of  observat ions.  The 
conventional  procedure  is  to  sort  the  output  samples  into  preselected  bins 
that  cover  the  range  of  interest  in  the  variable  y.  This  operation  of  sorting 
results  In  a  histogram,  which  can  be  integrated  to  form  the  sample  distribu- 
t  ion  f  nix't  ion. 

For  the  purpose  of  examining  the  upper  tail  of  the  distribution  function 
it  is  more  convenient  to  work  with  the  complement  to  P(y),  namely 

IX) 

C Hy )  ’  1  -  >’(y)  -  J 

y 


r(U  dt 


(l) 
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In  estimating  Q(y)  with  a  finite  set  of  observations,  we  can  define  the  oper¬ 
ation  applied  to  each  sample  as 


Dy(y)  -  1  .  y  >  Y 

-  0  .  y  <  Y  ,  (2) 


where  Y  is  a  preset  threshold.  In  practice,  this  operation  will  be  applied 
to  many  values  of  Y  simultaneously,  for  each  sample  of  y,  in  the  prior  step 
of  computing  the  histogram.  Ue  note  that  all  samples  receive  the  same  weight. 

The  result  of  applying  the  operation  in  (2)  to  each  sample  of  y  is  also 
a  statistical  process.  The  mean  value  of  (2)  is 

by(y)  •  j Dy (y)p(y)  dy 

-  /  P(y>  dy  »  Q(Y)  ,  (3) 

Y 


and  the  second  moment  Is 


d;<y>  « 


k 


(y)p(y)  dy 


oo 


■/ 


p(y)  dy  =  Q(Y) 


Tho  v.u  iiuiri'  is  j»ivon  I«y 


(4) 


var{  Dy  (y)  ]  =  Dy(y)  -  lDy'(y)]2 

=  Q(Y)  -  Q2(Y)  =  Q(Y)P(Y) 


(5) 


—  I  1i  III  Hi 


. 
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In  tin-  following  analysis  vc  will  bo  lnti-rcstcd  in  the  upper  tall  of  the  dis¬ 
tribution  where  F(Y)  ~  1,  so  we  can  essentially  assume  that 

var  [DY<y)l  =  Q(Y)  .  (6) 

In  order  to  estimate  Q(Y)  the  operation  in  (2)  will  be  repeated  for  N  samples 
of  y  and  the  estimate  will  be  formed  as 

A  1  " 

QO>  =  s  £  VV  •  (7) 

i  =  l 

where  {y^}  is  the  set  °f  N  observations  (statistical  samples),  which  we  assume 
to  be  independent.  It  follows  from  (3)  and  (7)  that  Q(Y)  ■  Q(Y),  which  means 
that  (7)  is  an  unbiased  estimate.  The  variance  of  (7)  is  given  by 

var  [Q(Y)]  -  ?  Q(Y)P(Y) 

=  iS,  Q(y)  •  (8) 

Tn  order  to  estimate  Q(Y)  with  big))  precision,  the  standard  deviation  of  the 
estimate  must  he  small  compared  with  the  mean  value.  In  other  words,  NQ(Y)»1. 
If,  for  example,  Q(v)  =  10  ^  then  N  must  be  at  least  as  large  as  10^  in  order 
to  achieve  any  precision  at  aLl  in  the  estimate.  This  requirement  for  an  im- 
pract really  large  number  of  samples  is  the  dilemma  faced  when  one  applies  con¬ 
ventional  sampling  techniques  to  the  estimate  of  low-probability  events. 

Importance  sampling  will  be  a  solution  to  this  dilemma,  but  before  we 
jump  to  that  subject  lot  us  expand  our  discussion  to  include  a  description  of 
the  process  itself.  An  sketched  In  Figure  1,  the  input  to  some  processor  will 
he  a  random  variable  x  with  a  known  probability  density  function  p(x).  The 
output  Is  y,  for  which  we  wish  to  estimate  Q(y),  the  complement  to  the  distri¬ 
bution  function,  with  a  finite  set  of  obsrjrvat  ions .  The  estimate  will  be 
formed  by  faking  the  average  value  of  the  observable  z,  which  is  the  operation 
by(y)  oppli.’d  to  y.  The  transfer  function  of  the  processor  will  be  designated 
by  y  K(x),  which  implies  that  a  given  value  of  x  is  mapped  into  a  unique  value 
of  y.  Since  each  x  lesults  in  a  paitlculac  value  of  z  we  can  write 


6 


Threshold,  Y 


input ,  x 
p(x) 


*  -  - 

rr«»cvssor 

output ,  y 

Kx) 

p(y) 

„ 

Sorting 
0  pi-mt  ion 
ny(y) 


observable,^ 


Figure  1.  Skrl  i'll  of  Opel  at  Ions  Per f 01  med  with 
Convent  ional  Sampling 
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*  *  DylK(x)| 


(9) 


2.2  IMPORTANCE  SAMPUHC  THEORY 

The  principle  of  Importance  sampling  is  to  distort  or  modify  the  input  ran¬ 
dom  process  (x)  in  order  to  make  the  original  low-probability  events  occur  more 
frequently.  This  action  will  be  compensated  by  weighting  the  event  by  a  factor 
that  is  a  function  of  only  the  particular  value  of  x  on  input  [1,2,3].  The 
functional  flow  is  sketched  in  Figure  2  where 

z„  -  Dv[F(x))w(x)  (10) 

10  I 


The  weight  is  designated  as  w(x)  and  the  modified  probability  density  function 

on  input  is  designated  as  p^(x).  In  order  for  z^  to  be  an  unbiased  estimate 

of  0(Y)  we  must  have  z  -  z,  the  latter  quantity  being  the  mean  value  of  (9), 
in 

so  that 


1 1  (x)  ]v(x)Pni(x)dx 


/' 


by  I F (x) 1 p(x)dx 


(11) 


thus 


w(x)  -  P(x)/P  (x)  .  (12) 

ID 

For  each  replication  of  the  experiment  a  specific  value  of  x  will  be  generated. 
With  the  use  of  (12)  we  can  then  compute  the  weight  on  the  basis  of  the  ratio 
of  known  input  probability  density  functions. 

So  far  we  have  been  working  with  a  one-dimensional  process.  Actually 
the  process  could  be  multi-dimensional.  We  will  assume  that  there  is  still  a 
single  output,  y,  so  that  we  can  write 

"  -  DylKfXj,  ....  xK)J  (13) 

and 

7m  *  Dy*F(x1>  •••»  XK>]U<*1 •  •••»  XK>  •  (14> 
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And  since  t  =■  z,  we  have 
m 


w(Xj  , 


*„) 


p(Xj ,  . . . ,  xK) 


(15) 


where  the  right  side  of  (15)  is  the  ratio  of  two  joint  probability  density 
functions.  If  the  input  samples  are  independent  and  belong  to  the  same  dis¬ 
tribution  function. 


K 

wCXj,  ....  xR)  =>  p<xv^pn<xit)  •  (16) 

2.3  SPECIFIC  EXAMPLES 

In  order  to  demonstrate  the  utility  and  power  of  importance  sampling, 
several  examples  will  he  given.  First,  we  will  examine  linear  processors 
where  the  output  stalistics  are  known,  and  then  we  will  analyze  some  non¬ 
linear  processors. 

Ex tmmlc  1:  Exponential _  Die tvibntion 

Let  x  be  exponentially  distributed  with  y  »  F(x)  -  x.  We  can  define 


p(x) 


(l/x)e 


-x/x 


x  >  0 


=  0  ,  x  <  0 


(17) 


where  x  is  the  mean  value  of  the  original  distribution.  Upon  integration 
we  obtain 

Q(x)  *  e‘x/*  (18) 

for  x  >  0.  For  importance  sampling  we  will  modify  t  lie  input  distribution  by 
(hanging  (increasing)  the  mean  value  as 


)0 


I>  (x)  -  2—  c“x/xm  ,  x  >  0 
is  x 


,  x  <  0 


The  weight  is  now  given  by 


w(x)  =  "  e-(1/x  '  1/xm)x 

P  (x)  X 


which  is  a  simple  calculation  that  will  be  performed  for  each  sample  on  input. 
For  the  case  when  y  -  F(x)  -  x  as  we  have  here,  the  second  moment  of  z 

a 

can  b<*  evaluated  easily.  For  a  single  observation 


"1  «  / Dy  <x)w2  (x)I’  (x)  dx 

mfi  m 


*(x)P(x)  tlx 


For  the  exponential  distribution  we  can  substitute  (17)  and  (20)  to  obtain 


-(2/x  -  1/x  )Y 
e  m 


x(2  -  x/x  ) 


Since  *  Q(Y)  -  c  ,  tlic  variance  for  a  single  observation  is  given  by 


var(-  | 


X|D  Y/x  .  1 

”  e  m-1 

x  (2  -  x/x  )  J 

m 


For  K  independent  observations  in  estimating,  the  distribution  function,  the 
vat  lance  will  be  reduced  by  a  factor  of  N  over  that  of  a  single  observation. 
Tine:  we  <  ..n  write  tire  following  ratio  for  N  independent  obsorval  ions : 
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var I “—] 

TD 


l 

N 


x(2  -  x/'x  ) 

ID 


Y/Xm 

TD 


(24) 


When  we  work  with  the  expression  in  (24)  It  Is  convenient  to  simplify 
it  slightly  by  anticipating  the  result  »  x.  Thus 


war  I  3 

m 


i 

N 


.  (25) 

» 


Now  wo  can  ask  the  question  as  to  what  the  optimum  value  of  x  is.  If  we  dif- 

_  m 

ferentiate  (25)  witli  respect  to  x  and  set  the  result  equal  to  zero  we  obtain 

ID 

x  -  Y  as  th«  solution.  At  tin*  optimum. 

to 


varh-.J 

‘  i,y 


e  Y 

2  :: 

x 


-1 


(26) 


For  illustration,  suppose  we  arc  int crested  in  estimating  Q(Y)  in  the  neighbor¬ 
hood  of  Q(Y)  «*  10  With  (18)  we  obtain  x  =  Y  =  I3.8x  (which  means  that  the 

m 

approximation  made  in  obtaining  (25)  was  indeed  valid)  and  (26)  reduces  to 


var  1 7.  J 
w 


L7.8/N 


(27) 


In  Section  2  with  con  writ  latt.il  sampling  techniques  this  ratio  was  10^/N.  We 
have  thus  reduced  the  number  ol'  samples  required  to  obtain  a  specified  preci¬ 
sion  by  a  factor  ol  over  '>0,000.  With  N  1000  samples  tin-  relative  standard 
deviation  in  the  estimate  of  Q(Y)  will  be  13.3%  with  importance  sampling,  an 
acceptable  error  for  practically  any  application. 

In  Figure  3  we  compare  the  experimental  distribution  function  with  N  *  1000 


n 


*  _  _ 

and  the  exact  distribution  function.  The  value  of  x^  “  13. 8x  was  used  to 
first  compute  the  histogram,  which  was  then  integrated  to  obtain  the  sample 
distribution  function.  While  the  procedure  was  optimised  for  Z(Y)  -  10-^, 
we  observe  that  relatively  Email  errors  exist  throughout  the  range  in  proba¬ 
bility  that  covers  five  orders  of  magnitude.  Thus  the  procedure  is  relatively 
insensitive  to  the  precise  value  of  x^  (which  could  also  be  established  by 
examining  (25)). 


Example  2:  Gaussian  Distribution 

Let  x  be  a  zero-mean  Gaussian  random  variable,  again  with  y 


where 


P(x) 


1 

/Tao 


F(x)  -  x. 


(28) 


and 


Q(Y) 


<x> 

/ 


P(x)  dx 


(29) 


For  importance  sampling  we  will  modify  p(x)  by  changing  (increasing)  o  so  that 


P  (x) 

ID 


] 

/2n  O 

ID 


-x2/2o 

e 


2 

ID 


The  weight  is  given  by 


P(x)  °m  -(I/O2  -  1/o2)x2/2 

"Mi')  “b  e 

ID 


(30) 


(31) 


The  inputs  to  the  process  were  exponentially  distributed  random  variables. 
Tf  u  is  a  uniformly  distributed  random  variable  (0,1),  then  x  -  -ln(u)  generates 
an  exponential  random  variable  with  unit  mean.  With  importance  sampling  the 
piecise  Eilat  istieal  properties  of  the  pseudo-random  numbers,  u,  are  no  longer 
critical  as  they  would  he  with  conventional  sampling  techniques. 
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From  (21)  wo  obtain  for  a  single  observation 


2 


a 


(32) 


where  Q(Y)  is  given  by  (29).  The  relative  variance  for  N  independent  obser¬ 
vation  is  given  by 


war Iz  ) 

m 


1 

N 


o  /o  q(y/£  -  a2 fa2  ) 

ia  x  _ n  / 

A~oW  «2<Y> 

m 


(33) 


In  anticipation  of  the  result  o  »  o  we  might  be  tempted  to  simplify  (33) 

by  setting  /2  -  (j  /o  “*  however,  the  resulting  approximation  would  not 

yield  an  optimum  solul Ion  ns  a  function  of  o  .  As  it  stands,  (33)  is  best 
'  m 

bn  nd  1  ed  nume  r  i  c«i  1 1  y  . 


We  can  iiiwl  nn  opt  iinum  value  of  u  that  minimizes  (33)  for  a  specific 

ttt 

value  of  Y.  lor  example,  if  Y  “  4 . 7o  Q(Y)  -  10  .  A  value  of  o  -  4.8o 

m 

was  found  to  be  the  optimum  for  this  value  of  Y,  although  (33)  was  nearly  flat 

over  the  interval  4.0  <  o  /o  <  5.5  so  that  we  can  conclude  that  O  *  Y  is 
m  ~  m 

essentially  the  optimum.  For  °m  “  ^  *  4.7o,  the  relative  variance  is 


var|.  ) 
m 

(7  )2  *  49/N  ,  (34) 

m 

which  means  that  N  1000  observations  would  result  in  a  relative  standard 

deviar ion  of  22%  in  the  estimate  of  Q(Y).  In  Figure  4  we  show  the  comparison 

between  the  experimental  (N  =  1000  and  o  =  4.7a)  and  exact  distribution  func- 

ra 

lions  for  the  Kmissinu- •  .ise.  While  o  was  chosen  to  optimize  the  procedure 
-  6  ra 

for  Q(Y)  ~  10  ,  wc  note  that  the  error  i s  •  relat i vely  small  over  most  of  the 

range  in  probability  plotted  In  the  figure.  We  see  also  that  the  error  is 
somewhat  larger  than  in  Figure  3  for  the  exponential  distribution,  as  we  could 
pr«’di*  t  from  (?./)  anti  (34). 
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Practically  every  situation  of  interest  in  the  analysis  of  radar  and  com¬ 
munication  systems  can  be  handled  by  either  the  exponential  or  Gaussian  random 
variable  as  the  input  to  a  process  (no  counter  examples  are  known  to  the 
author)-  For  example,  the  log-normal  distribution  is  a  simple  transformation 
of  Che  Gaussian,  and  the  Welbull  of  an  exponential.  In  fact,  one  can  create 
an  exponential  random  variable  from  the  Gaussian  and  vice  versa  as  we  will  show 
in  the  next  two  examples. 

Examvle  3: _ !A m-Squai-e  of  Two  Oauusian  Variates 

I.et  x^  and  x^  be  two  independent,  zero-mean  Gaussian  random  variables  at 
the  input  to  our  processor,  and  let 

y  =  F(x j .x^)  =  xj  +  x2  •  (35) 


If  we  select  (2  *  0.r)  then  y  will  be  an  exponential  random  variable  with 
y  «  1.  With  importance  sampling  the  weight  applied  to  each  outcome  will  be 
a  f unction  of  only  x(  and  x^-  From  (16),  (28),  and  (30)  we  can  write 


.  .  ,  2  -(1  -  l/2o?)(x?  1  x*) 

w(x,,x.,)  -  2o  e  m  l  2 

i  /  at 


;int!  with  O’O 


,  v  0  2  -a 

w(x. ,x.J  -  2 o  o 
12  m 


1/2 o\)y 


(37) 


We  have  applied  the  above  procedure  for  generating  the  exponential  ran¬ 
dom  variable  and  with  importance  sampling  we  have  computed  the  experimental 

2 

distribution  function.  We  chose  o  =2.63  (or  2o  =13.8)  and  N  *  1000.  The 

m  in 

experimental  result  is  essentially  identical  to  that  in  Figure  3,  at  least  in 
t  In*  .it  isf  i«  i  I  sense. 

i\t  if'iflt'  A:  C.  nc: ration  of  Gauiwian  Variate  Jrom_  Exponent  ini 

I  .el  x  be  an  exponent  Ini  random  variable  at  the  input  to  our  processor. 
Within  t he  processor  we  will  create 


y  ~  F (x)  »  /2x  cos  6 


(38) 


7 


* 

where  0  is  a  uniformly  distributed  random  variable  (0,2ir).  The  random 
vai table  y  will  he  Gaussian  distributed,  and  if  x  -  1  the  mean  value  of  y 
will  be  zero  and  the  variance  unity.  With  importance  sampling  the  weight 
applied  to  each  outcome  will  be  a  function  of  only  x  in  this  example,  since 
the  generation  of  0  is  performed  with in  the  processor  and  will  not  be  con¬ 
sidered  as  an  input  variable.  The  weight  is  given  by  (20)  for  *  ■  1. 

We  have  generated  Gaussian  random  samples  by  means  of  (38)  and  have  used 
importance  sampling  to  estimate  the  distribution  function  in  the  range  of 
10'^  <  q(y)  <  10 The  result  is  stat i st ically  consistent  with  that  of  the 

~  _  2 

direct  method  in  Figure  4.  A  value  of  x  =  22.1  (•  4.7  )  corresponds  to  op- 

-6  m 

timizing  the  procedure  at  Q(Y)  =  10 

With  this  example  we  no  longer  have  a  unique  mapping  of  x  into  y,  which 
ought  to  suggest  other  possibilities  in  simulation.  In  addition,  we  can  con¬ 
figure  t lie  processors  so  that  all  inputs  are  either  exponential  or  Gaussian 
random  variables  as  we  will  show  in  Kxample  9. 

Kxamplr  5:  Hum  of  FJxponcniial  Variate  a 

The  sum  oi  K  independent  random  variables  that  are  exponentially  distri¬ 
buted  is  a  chi-square  random  variable  with  2K  degrees  of  freedom.  We  can  write 

K 

y  -  F(Xj,  ....  xR)  =■  £  xk  ,  (39) 

k“l 


whore  (x^  is  1  set  of  K  exponential  random  variables  on  input.  With  im¬ 
portance  sampling  the  weight  will  be  a  function  of  the  specific  values  of  x^. 
Kror^  (16) ,  (17),  and  (19)  wo  can  write 


w(x  ( ,  ....  XK)  --  (x 


n/x)K  IXP  f-°/x  *  t/xJ  ]£  xk) 

L  k  =  l  J 


-  (X  /x)K  «T(1/X  -  1/Xm)  y 


(40) 


The  use  of  (38)  with  sin  G  substituted  for  cos  0  will  produce  a 
second,  independent  Gaussian  random  variable. 


lfl 


In  Figure  5  we  show  the  experimental  distribution  obtained  by  using  im¬ 
portance  sampling  on  (.39)  for  K  *  5  and  N  =  1000  replications.  The  mean  value 

of  the  input  distribution  was  modified  as  x  -  4.7x,  which  optimized  the  pro- 

1  —6 

cedure  at  Q(Y)  “  10  .  Note  that  the  error  is  extremely  small  throughout  the 

five  orders  of  magnitude  in  probability  shown  in  Figure  5. 

Example  6:  Sum  of  Exponential  Variates  with  Arbitral Weight c 

For  the  previous  example  all  input  random  variables  were  uniformly 
weighted  in  (39).  Let  us  generalize  the  situation  by  assuming  arbitrary 
weights  in  the  summation  as 

K 

y  '  K<xi . V  ”  £  Yk  •  (41) 

k=l 


where  (x^)  is  the*  set  ol  independent  exponential  random  variables  on  input 
with  a  mean  value  x.  With  importance  sampling  the  weight  is  a  function  of 
only  the  input,  so  the  weight  remains  the  same  as  (40) .  The  distribution 
function  ot  y  is  given  by 


<My) 


i 

k=i 


K-l 

[‘k 


II 

ii<k 


a.) 


-y>av 


(42) 


In  Figure  6  we  have  used  importance  sampling  to  estimate  Q(Y)  for  K=2 

and  N  1000.  The  ratio  of  weights  used  in  (41)  varied  from  .48/. 52  (resulting 

in  o**;eiit  in  I  ly  a  chi-square  distribution  with  4  degrees  of  freedom)  to  .05/. 95, 

with  tin  sum  of  weights  being  unity  in  all  cases.  For  each  of  the  cases,  the 

pi  ori'dm  e  was  optimized  in  t  lie  neighborhood  of  Q(Y)  ~  10  ^  (e.g.,  x  =  10. 7x 

m 

for  tin*  . 2r)/./5  case). 

With  example  we  ran  create  the  following  interesting  situation:  let 

K  '2  so  that  the  weigh!  used  in  imparlance  sampling  in  (40)  will  be  based  on 
two  Input  v.it  !ables;  however,  let  one  of  the  weights  a^  in  (4l)  be  zero  so 
that  I  In-  out  pul  y  will  bo  exponent  Lilly  distributed.  Fmin  one  replica!  ion  to 
the  ot  he  i  we  ran  even  i  tit  errh.inge  the  without  affecting  the  distribution 
fins'  ion  of  y.  Does  impoi  lance  sampling  work  in  this  case?  Indeed  it  does. 
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as  we  show  In  FI  guru  7  »  with  N  *■  1000  and  *  13.8.  The  error  is  not  as  snail 
as  it  was  in  Figure  3,  but  this  could  be  compensated  by  increasing  N. 

Exampl c  7:  Sun  of  tog -Normal  Variates 

The  distribution  of  the  sum  of  log-normal  random  variables  is  difficult 
to  compute  by  conventional  numerical  methods.  Let  us  begin  with  an  exponen¬ 
tially  distributed  random  variable  as  one  of  K  inputs  to  the  processor  in 
Figure  1,  with  x  =  1.  We  will  then  generate  a  Gaussian  random  variable  of 
unit  variance  as  in  Example  A, 

gk  "  COK  6  •  (43) 

where  0  is  a  uniformly  distributed  random  variable  (0,27!) .  The  log-normal 
random  variable  is  generated  as 

°LRk 

*k  -  «  .  (44) 

where  o  is  the  standard  deviation  of  the  log-variate.  The  median  value  of 
will  he  unity.  Tin*  output  of  the  processor  will  then  be 

K 

y  - 1  lk  *  <45> 

k-l 

for  which  we  will  use-  importance  sampling  to  estimate  its  distribution  func¬ 
tion.  The  weight  will  be  given  by  (40). 

Because  of  the  complicated  processor  above,  there  is  no  straightforward 
way  to  (house  which  modifies  the  input  distribution  for  Importance  sampling. 
Therefore,  we  will  try  some  arbitrary  values.  In  Figure  8  we  show  the  exper¬ 
imental  distribution  functions  for  eight  eases  where  we  have  varied  x  from 

m 

2  tlomie.li  30  ill  steps  Ilf  4.  In  each  case  K  r  2,  0^  ’  1.0,  and  N  ■  1000  ob- 
sci  vat  inns.  We  sec  r  ii.it  all  hot  one  of  the  cases  tend  to  he  clustered  about 

a  straight  line  in  the  tioin.it  plotted.  The  exception  is  the  x  =  2  case,  which 

m 
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ts  too  small  to  produce  values  of  y  frequently  enough  In  the  range  of  inter¬ 
est.  The  remaining  cases  establish  the  validity  of  importance  sampling  for 
this  example  since  each  estimate  is  unbiased;  the  spread  of  the  cluster  is 

a  measure  of  the  standard  deviation  of  the  estimate.  By  Inference,  x  -  1 

m 

(which  is  conventional  sampling)  with  a  sufficient  number  of  observations  would 
also  produce  a  result  that  would  fall  within  the  cluter  in  Figure  8.  Almost 
any  value  of  x^,  especially  in  the  interval  6  <  <  30,  could  be  used  with 

importance  sampling  to  produce  an  acceptable  result.  However,  since  x  s  22 
optimizes  the  procedure  tor  a  single  Gaussian  random  variable  at  Q(Y)  =  10  , 

there  is  probably  no  requirement  to  exceed  this  value  in  the  general  case. 

We  could  also  repeat  the  above  procedure  whenever  K  or  0^  were  changed. 

We  have  also  used  importance  sampling  on  (45)  by  beginning  with  indepen¬ 
dent  Gaussian  random  variables  as  the  input  to  the  processor.  The  above  con¬ 
clusions  were  unaffected  by  this  change,  however  (the  range  2  <  o  <8  proved 

~  m  - 

sat isi act  or y) . 


K.nuni'i  r  H:  !'um  of  Loo -Normal  Fhaaorc 

A  situation  that  is  common  in  the  analysis  of  radar  clutter  is  the  sum¬ 
mation  of  random  phasors  where  the  amplitude  of  each  is  log-normally  distri¬ 
buted.  We  will  generate  the  log-normal  random  amplitude  £.^  as  in  the  previous 
sect  ion  (but  now  refers  to  the  signal  amplitude)  and  then  generate  the  ran¬ 
dom  phasor  as 


v 


k 


(46) 


where  is  a  uniformly  distributed  random  variable  (0,  2ii)*  Finally,  we 
i  or  in 

y  it v* 

lui  which  VI'  will  mm-  iin|itii  t  ance  samp  1  i  up,  in  estimate  Us  distribution  function. 

In  I'ii’.mc  0  vc  bow  I  lie  expel  ini.nl  ;il  il  i  si  r  I  but  I  on  functions  for  three 
cases,  xm  18,  22,  and  2b,  where  K-2,  1>L  ^  0.5,  and  N  ^  1000  for  all  cases, 
nil-  cxpci  iiiicnt  al  curves  .-ip.atn  cluster  fairly  closely  throughout  the  entire 
rally,*-  in  |>r  eKibi  I  i  t  y  |>1ntted. 


(47) 
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Example  9:  Mixed  Eta  tie  lice 

111  Kx.  iiup  I  i  s  ■)  tlirour.il  8  wo  li.ivo  (super  Imposed  random  sample!,  from  tlie  same 
statistical  distribution.  But  in  radar  and  communication  systems,  signals  are 
often  combined  with  mixed  statistics.  For  example,  radar  clutter,  which  is 
often  assumed  to  be  log-normally  distributed,  would  be  combined  with  thermal 
noise,  which  is  Rayleigh-amplitude  distributed.  In  this  example  we  will  sim¬ 
ulate  this  case  and  use  importance  sampling  to  estimate  the  distribution  func¬ 
tion.  We  will  again  begin  with  two  exponentially  distributed  random  variables 
as  the  input  to  the  processor  in  Figure  1.  A  log-normal  phasor,  V^,  will  be 
generated  from  one  of  the  exponentially  distributed  inputs,  x^,  by  means  of 
(43),  (44),  and  (46).  The  Rayleigh-amplitude  phasor  will  be  generated  from 
the  second  exponentially  distributed  input,  x^,  as 

V2  =  /x^  (cos  4>  +  j  sin  G)  ,  (48) 

where  0  is  a  uniformly  distributed  random  variable  (0,  2it)  that  is  generated 
within  tin'  processor.  The  output  of  the  processor  will  be 

y  =  ivj  +  uV2|2  ,  (49) 


where  a  Is  a  factor  used  to  scale  one  process  with  the  other.  We  will  con¬ 
strain  x  =  1  for  both  inputs. 

In  Figure  10  we  show  the  application  of  importance  sampling  used  to  es¬ 
timate  the  distribution  function  of  y  for  a  =  1  and  N  »  10*  observations.  The 
four  values  of  -  14,  18,  22,  and  26  were  used  to  generate  the  four  experi¬ 
mental  cuives,  Note  that  they  are  tightly  clustered,  especially  for  the  higher 
probabilities  shown  In  the  figure.  We  have  also  shown  the  exact  distribution 
functions  for  the  components  that  make  up  the  sura  in  (49). 

Since  the  four  experimental  curves  in  Figure  10  are  so  closely  clustered 
throughout  five  oidcrs  of  magnitude  in  probability,  even  with  such  widely  dif¬ 
ferent  values  of  the  importance  sampling  scaling  parameter,  x  we  can  conclude 

m 

with  high  ptobability  that  the  true  distribution  also  lies  within  the  cluster. 
Thus  we  can  also  conclude  that  importance  sampling  works  when  the  process 


?8 


under  investigation  involves  a  combination  of  random  variables  with  differ¬ 
ent  statistics.  In  the  application  of  importance  sampling  one  probably  should 
create  all  random  variables  by  some  transformation  of  a  set  of  input  random 
variables  with  a  common  distribution. 

Exanip  U  10 CFAR 

A  simple  approach  in  radar  to  achieving  a  constant  false  alarm  rate  (CFAR) 
in  the  presence  of  nonstationary  noise  is  to  set  the  detection  threshold  on 
the  basis  of  the  average  noise  power  in  a  number  of  reference  samples  where 
each  of  these  samples  is  assumed  to  represent  noise  only.  Such  a  scheme  is 
denoted  as  cell  averaging  CFAR.  The  output  of  the  processor  might  be 

y=-Kx-  (so) 

X  . 

1 

where  x  is  the  sample  under  test  and  the  set  of  samples  {x.}  form  the  CFAR 
reference.  Note  that  if  we  modify  all  R+l  input  variables  equally  by  changing 
the  common  tin '.in  value ,  the  output  y  remains  unchanged  and  importance  sampling 
would  seem  to  he  inapplicable  in  this  case.  One  could  apply  importance  sam¬ 
pling  to  y  if  p (y )  were  known;  however,  this  would  be  a  severe  limit  to  the 
approach. 

In  order  to  find  a  method  that  works  on  the  input  variables  it  was  de¬ 
cided  to  treat  the  numerator  in  (50)  as  the  input  process;  the  denominator 
will  be  generated  internally  to  the  processor,  just  as  we  generated  0  in 
F.xainpl v  U.  The  reasoning  behind  this  choice  is  that  the  variance  of  the  de¬ 
nominator  is  much  lens  than  I  he  numerator,  and  the  low  probability  behavior 
of  the  numeiator  has  a  much  greater  impact  on  the  outcome  than  does  the  de¬ 
nominator.  However,  unless  R  is  reasonably  large  we  will  not  get  much  of  an 
improvement  with  i nrpo rt  mice  sampling  compared  to  the  conventional  approach. 

If  lire  random  variables  are  exponent  in] 1y  dist ribut ed *  the  following  number 
of  (disci  vat  ions  N  will  he  needed  lor  each  value  of  R  to  provide  the  same 
quality  of  result  with  importance  sampling  lQ(y)  ~  10 
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R 

oo 

200 

100 

50 

20 

N 

1000 

2300 

4700 

15,000 

170,000 

In  most  applications  of  cell  averaging  CFAR  in  radar  we  will  also  employ 
some  nom-olievent  integration  of  samples  prior  to  forming  the  CFAR  ratio.  In 
such  cases  w«*  can  write  the  output  as 


y 


R 


(51) 


Usually  then-  will  be  M  reference  ij.nnplcs  for  each  sample  in  the  numerator* 
so  we  can  write  R  =  KM.  In  a  typical  situation  we  might  have  K  =»  6  and  M  ■  10 
so  that  R  «  {>0,  which  permits  us  to  use  a  reasonably  small  number  of  replica¬ 
tions  to  determine  Q(Y)  with  importance  sampling  applied  to  only  the  numerator 
of  (51).  In  Figure  11  we  show  the  results  of  this  case  with  N  *  10,000  obser¬ 
vations;  all  random  variables  are  exponentially  distributed  and  x  =  5x  for 

ID 

the  K  »  6  input  variables.  Thus  importance  sampling  also  works  with  cell 
averaging  t'.FAR,  but  we  have  been  forced  to  redefine  the  input  to  the  processor. 

2.4  CONCLUSIONS 

Import  .nice  sampling  has  wide  application  in  the  simulation  of  signals  in 
radar  and  comnunJral ion  systems.  It  is  robust  and  efficient,  producing  re¬ 
liable  estimates  of  the  low-probability  tail  of  the  distribution  function  with 
typically  1000  replications  of  the  experiment.  It  has  been  shown  to  work  with 
multiple  input  processes  when  all  Inputs  belong  to  the  same  distribution  and 
are  distorted  equally.  For  a  process  involving  a  combination  of  several  types 
of  signals,  each  f tom  a  different  statistical  family,  one  can  still  apply  im¬ 
portance  sampling  by  redefining  the  procedure  in  which  the  signals  are  gen¬ 
erated.  In  effect,  a  new  processor  is  created  so  that  all  input  signals  will 
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belong  to  the  same  statistical  family  and  all  will  be  distorted  equally  in  in- 
portance  sampling.  In  S()iiie  cases,  such  as  with  cell  averaging  CFAR,  we  will 
even  have  to  redefine  the  input  to  tne  processor.  It  has  been  shown  that  one 
can  begin  with  either  Caussian  or  exponentially  distributed  signals  and  gen¬ 
erate  practically  any  combination  of  statistical  signals  used  in  the  simulation 
of  radar  and  communication  systems.  We  are  unaware  of  any  situation  of  practi¬ 
cal  interest  where  this  procedure  would  fail. 
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3.  MODULATION  OF  CLUTTER  WITH  A  SCANNING  ANTENNA 


There  are  two  basic  ways  to  scan  the  antenna  in  a  ground  based 
search  radar:  cither  continuously  or  step-wise.  In  the  latter  case 
we  can  assume  that  there  is  no  antenna  motion  relative  to  the  ground 
for  a  fixed  time  which  we  designate  as  the  on-target  time.  Received 
pulses  from  the  ground  clutter  will  be  partially  correlated  which  is 
due  entirely  to  the  internal  motion  of  the  clutter  itself  (the  wind- 
induced  motion  of  trees,  etc.).  If  the  antenna  moves  from  pulse  to 
pulse,  then  the  return  fron  any  point  on  the  ground  will  be  amplitude 
modulated  as  a  result  of  the  time-varying  antenna  gain  in  the  direction 
of  that  point.  The  correlation  properties  of  ground  clutter  will  then 
be  a l fee tori  by  the  scanning  modulation  in  addition  to  the  internal 
motion  of  the  duller.  In  many  cases,  the  scan  modulation  effect 
dominates  the  internal  motion  effect,  and  for  simulating  clutter  signals 
the  scan  modulation  severely  complicates  an  otherwise  straightforward 
impl emen t  a t ion . 

For  the  purjK>ses  of  this  discussion  we  will  distinguish  between 
two  types  ol  radar  signal  processing:  continuous  and  batch.  With 
continuous  processing  the  radar  generally  utilizes  all  past  data  in  a 
recursive  filter  to  process  each  pulse,  and  there  is  usually  one  output 
signal  for  each  pulse.  Continuous  processing  is  often  Implemented  in 
continuously  scanning  systems,  at  least  for  older  radars.  With  batch 
piotessiug  the  radar  collects  a  sequence  of  pulses  before  any  processing 
is  implemented,  and  the  processing  from  one  batch  to  the  next  is  independent. 
Hatch  processing  is  implemented  on  all  step  scan  systems  (as  known  to 
the  author)  and  some  of  the  newer  continuous  scan  systems,  especially  the 
3D  search  radars  utilizing  an  electronic  scan  in  elevation.  The  FFT 
processor  is  one  example  of  a  hatch  processor. 

The  distinction  between  continuous  and  batch  processing  has  a  dra¬ 
matic  effect  on  the  techniques  used  to  simulated  clutter  signals,  so  much 
sc*  that  the  techniques  and  computer  algorithms  used  for  one  system  will 
bear  little  similarity  to  those  used  for  the  other.  In  effect,  there 
will  be  two  distinctly  different  computer  programs  required,  to  simu¬ 
late  both  « uses.  The  reason  is  that  the  continuous  processing 
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system  demands  almost  a  pulse-by-pulse  approach  to  computing  the  inter¬ 
action  of  the  antenna  with  the  clutter  geometry,  while  the  batch  pro¬ 
cessing  system  usually  needs  to  be  updated  only  once  for  each  batch  of 
pulses.  The  latter  case  is  simpler  and  cheaper  to  simulate,  and  few 
questionable  approximations  will  be  required  out  of  necessity.  For  this 
reason  we  will  concentrate  our  discussion  on  the  continuous  processing 
system.  We  will  assume  in  the  following  analysis  that  the  radar  has  2-D 
resolution  and  the  antenna  has  a  fan  beam  in  elevation. 

3.1  FORMULATION  OF  THE  CLUTTER  SIGNAL 

Let  us  designate  the  complex  signal  associated  with  the  k*'*1  scatterer 
in  a  range  resolution  cell  as  V^Ct),  which  we  have  made  time-varying  in 
order  to  account  for  internal  motion  of  the  clutter.  In  general,  |v^(t)|^ 
will  have  the  dimensions  of  radar  cross  section  (RCS),  or  some  power 
scaled  version  of  it  according  to  the  radar  range  equation.  In  addition, 
let  g(0)  be  the  one-way  voltage  gain  of  the  antenna  in  the  azimuth,  8, 
direction.  The  composite  signal  (complex  voltage)  received  by  the  radar 
when  the  antenna  scans  continuously  in  azimuth  will  be  proportional  to 


V(t)  =  Svk(t)g2(eo+0t-ek)  (52) 

k  " 

where  0^  is  some  reference  angle,  0  is  the  scan  rate,  6^  is  the  azimuth 
of  the  k1*1  seat  I  erer,  and  the  .summation  is  performed  over  all  scatterers 
within  a  range  resolution  cell.  Differences  in  range  among  scatterers 
that  would  etfecl  the  phase  of  superposition  in  (52)  are  implicitly  included 
In  the  definition  of 

As  written  (52)  Is  nontrivial  to  implement  because  vk(t)  is  a 
two  dimensional  process :  spatial  and  temporal.  Without  any  simplifica¬ 
tion  (which  can  be  done  only  through  approximations)  we  must  evaluate 
(52) lor  each  pulse,  and  at  the  very  least  we  roust  have  one  scatterer 
(in  a  range  ling)  for  each  pulse  throughout  a  360”  scan.  For  a  PRF 
of  A 00  llz  and  a  10  sec  scan  wc  will  have  4000  total  pulses  and  at 
least  4000  .scatterers.  The  process  V^(t)  is  thus  described  by  a 
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4000  x  4000  matrix  (neglecting  the  "edge  effects"  at  the  beginning 
and  end  of  the  scan).  Clearly,  some  kind  of  approximation  is  in  order 
at  this  point  to  reduce  the  size  of  the  problem. 

One  approach  would  be  to  assume  that  g(9)  is  negligibly  small  for 
1 0 1  larger  than  some  value  (where  0*0  corresponds  to  the  main lobe 
axis.)  We  could  simulate  only  the  mainlobe  or  the  near  sidelobes  of 
the  antenna,  and  assume  that  there  is  no  contribution  from  sidelobes 
furts.fr  out  than  some  point.  The  problem  is  that  in  some  cases  (in  older 
systems)  just  simulating  the  mainlobe  might  be  adequate  but  in  other  cases 
(in  newer,  high-performance  systems)  we  might  have  to  include  many  sidelobes. 
Just  where  we  should  slop  simulating  sidelobes  is  not  readily  apparent 
since  it  depends  on  the  system  being  simulated.  Few  ground  rules  can  be 
given  either  because  very  little  analysis  has  been  performed  on  this 
subject . 

Another  approach  to  simplifying  (52)  is  to  focus  our  attention  on  V^(t). 
In  general,  this  process  can  be  assumed  to  be  random  (for  the  objective 
of  determining  detection  per  I onnance) ,  with  the  statistics  of  one  scatteror 
being  independent  fioiu  .my  other.  The  random  process  will  be  spatially 
uncorrelated.  It  is  fairly  common  in  the  literature  to  break  up  the 
scattering,  properties  of  ground  clutter  into  two  components:  an  ac 
component  and  a  dc  component.  The  ac  component,  the  time-varying  one,  is 
the  result  of  t'..r  motion  of  trees,  brush,  grass,  etc.,  and  is  generally 
made  up  of  many  individual  scattercrs  (e.g. ,  the  leaves  on  a  tree) 
whotr  no  single  one  dominate.'**  Such  a  process  is  easily  justified  to 
be  C.aussian,  at  least  within  a  relatively  small  area  on  the  ground. 

On  the  other  hand,  the  dc  component  is  associated  with  rigid  objects, 

•  uch  as  ha i e  ground,  tree  tnniks.  man  made  structures,  etc.,  and  non- 
<;.nr.*;  ian  statistics  usually  prevail  in  this  case.  In  general,  the  dc 
component  is  stronger  (higher  power)  than  the  ac  component. 

With  the  above  distinction  hetwe«>n  the  two  component  s  let  us  define 

V>  *  \  1  V«) 

wli.1,-  Vk  is  t  lie  lie  . . .  of  I  lie  k,h  acnttetci  .Hid  Vfc(t)  is  the  ac 

io".|>on.iit  .  It  w.-  ;.ul.  I  1 1  ute  (7)  into  (1)  we  obtain 


(53) 
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vu)  -  HvlR2<« +et-o.)  +  ]Cv.<t)g2<e  +ot-o.)  (54) 

k  K  °  *  K  ok  j. 

At  iirst  glance  it  appears  that  we  have  complicated  the  situation  because 
we  now  have  two  summations  insted  of  one.  However,  the  first  one  is 
straightforward;  it  is  a  circular  convolution,  something  that  is  easily 
implemented.  For  the  second  summation,  which  has  all  the  inherent  complexity 
of  (52)  we  will  truncate  the  antenna  pattern  as  described  previously,  but 
now  the  effect  of  the  truncation  will  be  much  less  significant  since  V^(t) 
will  usually  be  much  weaker  than  V^(t). 

Equal ion  (54),  or  at  least  the  second  term  of  it,  is  not  particularly 
economical  to  implement  in  terms  of  computer  time,  especially  when 
computing  fnlse-alarm  performance.  For  this  reason  we  will  pursue 
further  approximations. 

3.2  FURTHER  APPROXIMATIONS 

In  general  (52)  and  (54)  desci ibe  a  nonstationary  process  because  clutter 
samples  in  one  azimuth  sector  can  be  much  stronger  than  in  another.  This 
is  another  way  of  saying  that  clutter  is  spatially  nonhoraogeneous.  If  this 
property  of  clutter  applies  lo  both  terms  in  (54)  then  there  is  little  more 
we  can  do  •«'  simplify  the  expression.  However,  if  the  ac  component  of 
clutter  is  homogeneous  (more  speculat ion  at  this  time)  we  can  treat  the 
second  term  in  (V*)  as  .a  stationary  process.  With  this  assumption  we 
can  apply  Fourier  analysis  to  determine  the  spectrum  of  this  term  which 
we  rewrite,  as 

v.ic(t)  -  L.  vk(t)g2(eo+0t-0k)  (55) 

In  determining  the  second  moment,  let  us  write 

V.-,c(t.)Vl(t2)  •  pk(tl)Vt2)S2(V»l-Ve2(V't-9l)  (56) 

Wln*(i  we  take  the  «*usemhle  average  of  (56)  the  cross  terms  vanish  as 

WV‘2)  “  0  •  y  *  1 
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Vac(ll>VI>2>  *  ?  B2(Bo+etr9k)g2(9o46t2-6k) 


autocorrelation  function 
of  composite  process 


autocorrelation  function  of 
clutter  internal  motion 


The  autocorrelation  functions  above  will  be  assumed  to  be  nonstationary 
(even  though  we  will  later  assume  that  the  autocorrelation  function  of  the 
clutter  internal  motion  is  stationary).  Thus  we  will  define 

“(,l-'2>  =  VV'frV  <58> 

Rac(Vt2)  "  Vac<tl)Vac(t2>  (59) 


hb(t)  -  g  (6  +  Ot  -  0.) 

K  OR 


Note  that  we  are  assuming  that  (58)  is  independent  of  the  subscript  k  so 
that  the  clutter  process  Is  spatially  homogeneous.  With  this  notation,  (57) 


Rae<V2>  "  I  R<V‘2>  hk  (tl>  \  (t2> 


At  this  point  we  will  introduce  the  two-d imens ionnl  Fourier  transform 
as  nf.ed  in  Refer  ome  7,  so  that 

“  fj  c~j2l'(fltl  "  f2t2)  dt1dt2  ,  (6 1 


*8 


for  the  moment.  Then  it  can  be  shown  that 

XkU.O  -  VS  e-”2f2/2a  e-aT2/2  e-^fl  . 

Now  if  the  nntenna  pattern  is  defined  by 

g2(0)  -  e-<“0/63dBi2  . 

where  a  ~  1.6651  and  0^^  Is  the  one-way  half-power  hrnmvidth,  we  note  that 
from  (60)  and  (68) 

*  *  <a6/°-,dB)2 

and 

t  -  t  +  (0  -  0  )/0 

o  k 

Thus  (69)  is  Ivon  by 

x  o.f)  ,'(",0»UH/“")2e  '('■“l/D1dB>2f^,(Wi-1,2)  •  (73) 

k  *  2  up 

Note  that  in  (67)  the  subscript  k  is  only  associated  with  the  ambiguity 
fiiin  tinn,  and  in  (7  1)  the  suhueript  k  only  appears  in  the  last  term.  Thus 
it  we  sum  (73)  over  k  we  must  ova  Info  the  following  expression 

l  <J2vi  <<V0k)/®~T/2)  *  (74) 

k 

lli'W’Vi-i  ,  .»•:  tin*  inen  menlH  become  small  so  that  the  summation  approaches  an 
i  nt  ey,r  at  ion,  (74)  will  have  a  significant  value  only  for  f  ~  0.  As  the 
limits  on  the  Integration  get  larger  and  larger,  (74)  becomes  more  like  a 


(71) 

(72) 


(69) 


(70) 
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6-funrl  Ion,  6(f).  For  our  purposes  we  will  assume  Hint  it  is  a  6-function 
bemuse  we  will  never  employ  long  enough  integration  lines  Lo  be  nble  to  re¬ 
solve  the  difference  in  shape  of  (74)  from  a  6-function.  We  must  also  in¬ 
clude  a  scale  factor  in  the  evaluation  of  (74)  as 


e  eJ2„f((0o-0k)/0-,/2)  0  (75) 

~  AU 

k 

where  AO  is  the  azimuth  spacing  of  the  clutter  samples  around  the  range  ring 
(or  the  average  spacing  if  the  samples  are  nonuniformly  spaced).  With  this 
assumption  the  suumation  of  (71)  over  k  becomes 


Fiiitluimoro.  (67)  reduces  to 


(76) 


r  (f 

ac 


j.f2)  «  <5(1  |-f 2>  e~^2llf2T  dl  , 


where 


X(1)  -  J\ 


°3dB  -^(uOl/B  )2 
e  3dB 

uAO 


on 

(78) 


Note  that  (/7)  ir.  just  the  aut  ocorrel  at  ion  function  of  the  beam  scan  modu¬ 
lation  function  in  (60)  with  the  (.mission  beam  assumption.  Furthermore, 
(//)  is  ?.  i  np.nl  ar  alonp,  I  lie  fj  =  f,,  axis,  wliicli  means  that  t  lie  resultant 
process  is  slal  ionary  l/].  if  we  designate  ^ac^O  as  the  power  spectrum 
of  this  stationary  process,  then  we  can  write 

sacu)  - /  R(i)xd)  c"J2llfT  di 

Now  let  ns  define  the  following  Fourier  transform  pairs 
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SlB(f)  —  R(«)  (80) 

G(f)  -*•  X(r)  (81) 

so  that  the  spectrum  of  the  composite  process  in  (79)  can  be  written  as 


s  tf)  -  s.  ★c(f) 

ac  lm 


(82) 


where  S^m(f)  is  the  spectrum  of  the  clutter  Internal  motion.  G(f)  is  the 
scan  modulation  spectrum,  and  the  star  designates  a  convolution. 

l.et  us  continue  with  the  assumption  of  a  Gaussian  beam.  We  will  define 


f 

g 


a26 


/2n  0 


3dB 


(83) 


so  that  when  we  Lake  the  Fourier  transform  of  (77)  we  obtain 


«(0 


/v  °3dB 
2  "(iAO 


u  -(af/f  )2 

/a  i  e  8 

g 


(84) 


normalized  to  unit  power 


Note  that  it.  the  two-sided  ha If-power  spectra]  width.  Let  us  also  assume 
that  R(t)  is  a  ( Russian -shaped  autocorrelation  function 


R('0 


Pta_!  -  0>  t ) 2 

2tr  e 


(85) 


Note  that  this  autocorrelation  function  applies  to  a  single  clutter  sample; 
therefore  we  have  chosen  as  the  average  power  associated  with  all  clutter 
within  a  range  ring  throughout  2ti  radians  in  azimuth.  The  Fourier  transform 
of  (8!i)  is  also  C.tussian  shaped  ns 


S . 

i  m 


<f> 


iy\o 

2  ii 


-(‘xf/f.  y 


/h!  . 


normalized  to  unit  power 


(86) 
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where 

b  -  itf  /a 

IB 


and  f.  is  the  two-sided  half-power  width  of  the  internal  notion  spectrin, 
in 

The  convolution  of  (84)  and  (86»)  is 


S 

ac 


<f) 


/«  (l3dB 
/2  •  2* 


u  -(af/f  )2 
-  e  ac 

/a  f 


normalized  to  unit  power 


(87) 


I*  2  -way  half -power  bcarawidth  •  2ti 
1 ,064 

Iota!  power  ill  range  ring 


where 

lJ  ■=  i]  +  f2  (88) 

ac  iro  g 

We  can  neglect  the  scan  modulation  when  f.  >>  f  *  or  from  (83)  when 

im  g 


. 624  0/0 


3dU 


(89) 


We  can  replace  the  •  sign  with  a  >  4.0  without  any  noticeable  effect  on  per¬ 
formance  • 


42 


4.  ci:ni:katin<;  currrKR  skqukncks  for  a  ground-based  radar 

Ln  Reference  H  several  general  techniques  are  described  for  generating 
clutter  sequences.  They  are  all  based  on  properties  of  the  discrete  Fourier 
transform  and  the  fact  that  samples  in  the  frequency  domain  will  be  statis¬ 
tically  Independent  of  each  other.  We  will  generate  a  set  of  independent 
random  phasors  (X(nAf))  at  uniformly  spaced  increments.  Af,  in  the  frequency 
domain,  such  that 


|X(nAf)|2  =  AiS(nAf)  (90) 

where  S(f)  is  the  desired  power  spectral  density  and  the  bar  on  the  left  side 
of  (90)  designates  an  ensembel  average.  The  amplitude  distribution  of  the 
phasors  {X(nAf)J  need  not  be  Rayleigh  18 J.  The  correlated  time  sequence  is 
obtained  by  taking  the  discrete  Fourier  transform  as 


x(kAt)  X(nAf)  cJ2l,knAfAt  (91) 

n 

Usually  the  sample  spacing  in  the  time  domain.  At,  will  be  given  (e.g.,  the 
pulse  repetition  period  in  a  pulsed  radar),  but  Af  is  under  our  control.  Be¬ 
fore  we  define  how  Af  should  be  cliosen  we  note  that  the  time  sequence  repeats 
with  a  period  given  by 

Tr  «  1 /Af  (92) 

hot  us  also  del  ine  the  two  sided  half-power  width  of  the  spectrum  ns  f-jdB* 

The  correlation  distance  is  then  approximately 


T 

c 


1/f 


3dB 


(93) 


Since  the  tune  sequence  repeals  with  a  period  T^.,  we  can  utilize  only  a 
p«*i  l  ion  of  the  pel  iod  without  having  the  beginning  ot  the  desired  sequence 
being  correlated  with  the  end.  We  should  therefore  choose 

1  T  '  T  T  (94) 

r  “  c 


where  T  is  the  duration  of  the  desired  sequence.  If  we  make  use  of  (92)  and 
(93)  we  can  rewrite  (94)  as 


Af  < 


3d  It 


-  1  + 


f  T 
3dB* 


(95) 


There  is  one  other  constraint  on  Af :  We  must  be  able  to  resolve  the  8pectrua. 

From  Reference  8  (p.  120)  we  must  have  Af  i  0.63  f  for  a  Gauss  fan-shaped 

Jan 

spectrum.  The  constant  will  be  slightly  different  for  other  spectral  shapes* 
but  it  will  probably  not  be  loss  than  0.63  for  ground-based  radar  clutter 
spectra.  Therefore  wo  can  write 


Af  <  0.63  f .]  | u  •  min 
-  3dB 


1.60 


1  +  f3dBT 


(96) 


In  general,  will  be  less  than  0.60  for  a  ground-based  surveillance 

radar,  so  the  simple  constraint  Af  <  0.63  will  apply  in  most  cases;  how¬ 

ever,  (96)  will  always  be  applicable. 

There  are  two  general  ways  to  implement  (91):  the  fast  Fourier  transform 
(FFT)  and  brute-force  approaches.  We  will  now  discuss  each. 


4.1  KFT  APPROACH 

l.et  us  define  the  repetition  frequency  as 


I  *  1/At  (97) 

Next,  let  us  divide  this  interval  into  equal  increments  so  that 

Af  *  f  /N  (98) 

r  r 

From  (97)  we  note  that 

Nr  1 /Af  At  (99) 

so  that  (91)  can  he  written  as 

Nr  *  j2nkn/N 

x(kAt)  -  }2  X(nAf)  «'  r  O00) 

n*U 
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which  is  the  conventional  definition  o!  the  discrete  Fourier  transform  that 
can  be  implemented  as  an  FKT. 

In  order  to  determine  how  large  should  he  let  us  utilize  (98)  and  the 

inequality  Af  <  0.63  f  ,  and  write 
JQD 


Nr  :  1.60  fr/f3dB  «01) 

with  tin-  understanding  that  might  have  to  be  even  larger  if  (95)  were  to 
apply.  Now  let  us  work  with  some  numbers.  A  general  rule  of  thumb  is  that 
the  spectral  width  ol  ground  clutter  for  a  ground-based  radar  with  a  non- 
scanning  antenna  is  about  'iZ  of  the  maximum  wind  Doppler  1$).  Thus  we  can 
write 

f  .  lu  -  .06  V  /A 
IdB  w 

win-re  V  is  the.  wind  velocity  and  A  is  the  wavelength.  For  V  *  10  m/sec  and 
w  J  w 

-  *  .  1/  m  (b-baud)  we  have  l ..  ,  -  5  II/..  If  the  PRF  is  1  kHz  (f  *  1000  Hz) 

3d  d  r 

then  (101)  becomes  N  j  320,  which  is  a  large  number,  especially  considering 
that  we  might  utilize  at  most  only  about  16  samples  in  the  time  domain, 
i  he  d  i  sadv.iul  ag<  ol  the  KFT  approach,  at  least  as  defined  so  far,  is  that 
most  of  the  spectral  samples  will  be  zero.  Kvon  though  the  FKT  fs  efficient, 
it  must  still  implement  the  multiplies  by  zero.  The  next  approach  avoids 
ibis  shortcoming. 

‘.2  HKim;-  FOKOF  AITKOACI1 

The  brute  force  approach  is  a  direct  imp  J  ement.it  ion  of  (91).  We  define 

l  lie  samples  in  the  f  liqueur y  domain  only  over  a  limited  region  of  the  power 

■  poettal  density.  II  we  utilize  Af  0.63  f  ,  ,  then  as  few  as  5  or  7  samples 

Job 

•ill  be  sufficient  to  define  the  spectral  process  (as  determined  in  Refer¬ 
ee  c  H  for  Caussi.m  shaped  spectra;  the  rule  might  be  different  for  other 
-.port  in  I  shapes).  Moicovei  (91)  is  implemented  only  for  the  desired  number 
of  time  samples.  If  and  N  are  the  number  of  frequency  and  time  domain  sam¬ 
ple;.,  respectively,  then  the  computation  time  will  be  proportional  to  the  pro¬ 
duct  NNj  ,  instead  nl  N^log^N^  for  the  FIT  approach.  But  since  N  <<  and 
N j  ■  will  usually  be  true  for  ground-based  radars,  the  brute  force  approach 
will  usually  he  faster  than  the  FF'f  approach  for  generating  clutter  spectra. 
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4.3  FFT  APPROACH  WITH  INTERPOLATION 

If  '*•  f  jjj,*  as  it  is  with  most  ground-based  radars,  then  many  consecu¬ 
tive  time-domain  samples  will  be  correlated.  We  can  increase  the  sample 
spacing  At  in  (91)  or  (100)  to  reduce  the  computation  time  for  the  Fourier 
transform,  and  then  utilize  interpolation  to  obtain  samples  at  the  desired 
rate.  Let  us  define  the  reduced  sample  spacing  in  the  time  domain  as  At', 

.iiul 

h  =  At /At'  (102) 

Similarly,  we  will  define  c  1/At',  so  that 

h  =  f  '/f  (103) 

r  r 

l-ur ihermort* ,  we  will  deline  the  number  of  samples  in  one  repetition  interval 
in  the  f requeue v  domain  as 


N  '  *  f  '/At  (104) 

r  r 

.mil  as  a  result 

N  =  N  '/h  (105) 

r  r 

Interpolation  in  the  time  domain  causes  spurious  responses  in  the  fre¬ 
quency  domain  [fcj.  1 f  we  wish  to  hold  these  spurious  responses  to  50  dB  be¬ 
low  the  desired  clutter  power,  we  roust  choose  f^*  >  10  if  linear  inter¬ 

polation  is  used  (8,  l-xj.  8.61).  Now  we  can  use  (103)  and  write  h  >  10  ^^dE^r' 
liowcvcr,  if  this  ratio  is  greater  than  unity  we  will  not  implement  interpola¬ 
tion.  Thus  we  can  write  an  equality  that  will  provide  -50  dB  spurious  responses 
for  linear  inlet polation  as 


h  =  minUO  f.WB/fr.  1) 

ff  wo  use  Af  <  0.63  f_,„  in  (104)  we  can  also  write 
3d  IS 


N  '  >  1.60  f  ' 
r  -  r 


/f 


3dB 


(106) 


007) 


which  is  equivalent  to  (101)  applied  to  the  reduced  sampling  rate.  If  we  com¬ 
bine  (107)  with  (103)  we  obtain 
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N 


r 


1 .60  li 


f  H 


3UB 


008) 


.aid  with  (1 06)  we  have 


N  '  >  16  .if  f_,_/f  <  0.1  (109a) 

r  -  3dB  r 


■*  1.6  f  /f_.  ,  otherwise  .  (109b) 

r  jdc 

Note  that  (109a)  provides  the  largest  lower  bound  on  N^* ,  so  that  we  can  write 
lor  all  cases  '  2  16.  This  is  a  factor  of  20  less  than  in  the  example 
in  Section  4.1. 

Let  us  define  N  =  T/At.  Then  from  (93),  (94),  and  (97)  we  can  write 


N  >  N  +  f  /f,1K 
r  “  r  3dB 


(110) 


With  (105)  and  h  11  10  f_.„/f  ,  We  obtain  another  lower  limit  on  N  1  as 
>db  r  r 


N  ’  -  (N  +  10)h 


(HI) 


••iu.illy  for  a  ground  based  surveillance  radar  this  limit  will  be  lower  than 
{ 1 09.i )  . 

In  Table  1  we  list  the  perf  or  man  rte  that  can  be  obtained  for  various  param¬ 
eter  choices.  The  use  of  the  table  will  begin  with  a  specification  of  the  de~ 
.ired  spurious  spectral  response  level  and  accuracy.  For  efficient  computation 
Hie  lowest  value*  of  '  should  be  used,  although  some  consideration  should  be 

'iven  to  the  use  of  =  16, 32, 64, etc. ,  because  these  values  are  especially 

•  ficieni  KIT  si/.«-s.  We  will  be  given  N,  the  number  of  time  samples  being 
utilized,  and  the  ratio  We  will  then  compute 


point  ion.  Next  we 


h  = 


.  given  in  the 
will  (heck  the 


(fr*/fW  f  tdh 


table.  If  h  1 
constraint  on  N 


/f  012) 

we  will  not  implement  inter¬ 
in  (111),  which  we  rewrite  as 


N  <  N 

r 


7"  - 


(UV 


AH 


If  il  is  sat  is!  ic*d  then  we  r.ui  use  that  particular  option;  otherwise,  we  will 
have  to  use  a  larger  value  of  '  .  let  us  now  consider  the  example  *  -01 

and  N  =  32*  and  assume  that  Opt  ion  A  would  be  suitable  if  it  satisfies  the  con¬ 
straint  on  N.  We  obtain  h  ■  0. 10  and  the  * onsl raint  N  £  ISO,  which  works.  On 
the  other  hand,  suppose  that  f  /I ^  *  .OS  and  the  performance  of  Option  E 
were  the  minimum  atceptable.  N«*w  we  vtain  h  ■  0.80  and  N  *_  64,  which  works 
(Options  K  and  <;  would  also  work).  We  shuuhi  comment  that  when  h  is  close  to 
unity,  as  in  the  last  example,  intn  |h*  l.it  ion  may  not  offer  any  computational 
advantage  over  the  straight  KKl  appro-n  it.  We  could  have  chosen  Af  =*  f^/64  * 

.3125  f.,.,  ,  and  wo  eould  have  utilized  N  *  64  -  20  *  44  time  samples  without 
MB 

Interpolation,  which  would  have  i  mulled  in  a  slightly  faster  computation.  In¬ 
terpolation  pays  off  when  h  <  O.S,  especially  so  for  very  small  values  of  h. 

When  O.S  <  l»  <  1.0  the  issue  is  not  so  decisive  as  many  factors  come  into  play. 

We  can  compare  the  compiit  at  ion  times  tor  the  brute-force  and  FFT  (with  in¬ 

terpolation)  approaches.  For  the  former  the  computation  time  is  proportional 
lo  N  •  Nj- ,  where  Nj  is  the  number  of  samples  defining  the  power  spectral  don- 

•  It*.  For  the  1-FT  the  computation  time  is  proportional  to  N^'  log^  •  On 

•  nr  paiticular  computer  tin  hiule-force  approach  will  be  faster  when  N  sat is- 
1  it  s  tin*  i  ol  lowing: 

N  ’  N  ri  N  =  7 
r  t  f 

lh  N  i  !>  N  *  4 

y?  N  i  14  N  £  10 

■  •  >  i  ?  lf»  the  const  mint  s  of  N  •  A  or  N  <  5  will  not  usually  be  of  practi- 

al  interest  in  i.ul.ir  simulat  ions;  even  if  wo  were  interested  in  .simulating 
o  lew  samples  it  would  he  still  relatively  efficient  ( o  utilize  the  FFT  np- 

•  roach  with  interpolation.  For  N^'  -  32  there  might  be  a  few  situations  in 
vhich  the  brute-force  approach  would  offer  a  computational  advantage,  but  the 

dvantage  will  never  be  great. 
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5.  SIMULATION  PROGRAM 

A  very  convenient,  flexible*  and  computationally  efficient  Fortran  program 
has  been  developed  to  simulate  detection  performance  for  a  ground-based  radar 
against  airborne  targets.  The  quadrature  components  of  the  video  signal  are 
simulated  in  the  neighborhood  of  the  target*  and  the  received  signal  is  pro¬ 
cessed  in  the  receiver  in  the  same  manner  and  sequence  as  it  would  be  in  an 
actual  radar.  The  scattering  environment  is  described  in  a  statistical  manner, 
as  well  as  the  fine-scale  location  of  the  target  within  the  azimuth  beam*  range 
gate,  and  Doppler  filter.  By  repeating  the  run  with  independent  random  inputs, 
the  detection  statistics  are  accumulated  in  the  form  of  cumulative  probability 
versus  threshold  setting.  if  a  target  is  absent*  then  the  output  is  probability 
of  false  alarm  versus  threshold  setting.  Importance  sampling  techniques  can  be 
used  to  increase  the  efficiency  of  this  calculation. 

Practically  any  situation  of  interest  can  be  simulated  with  parameters  speci¬ 
fied  on  input.  Listings  of  the  program  as  given  in  Appendix  A  include  a  descrip¬ 
tion  of  all  input  parameters,  as  well  as  the  procedures,  assumptions,  and  limitations 
within  each  simulation  step.  The  main  program  (MAIN)  acts  as  a  driver,  in  which 
calls  are  made  to  subroutines  that  generate  specific  signals  such  as 

•  the  target  (TARGET) 

•  clutter  (CLIJTTR) 

•  receiver  noise  (NOiSK) 

in  addition  to  a  subroutine  that  processes  the  received  signal  (PROCES)  and  ac¬ 
cumulates  detection  performance  statistics  (DSTINL,DSTPNT  ,1)STFUN) .  The  simula¬ 
tion  procedure  can  be  easily  modified.  For  example,  as  written,  the  main  program 
permits  two  types  of  clutter  to  be  generated,  ground  plus  rain  or  ground  plus 
chaff;  all  three  types  can  he  generated  with  three  calls  to  CLUTTR  (which  also 
requires  additional  parameters  to  be  defined  on  input). 

5. 1  INPUT  PARAMETERS 

All  input  parameters  are  capable  of  being  defined  in  DATA  statements  in  the 
main  program  or  via  NAMKI.IST  /VALUES/  for  each  simulation  run.  We  will  discuss 
these  p.u  amet  iTS  as  they  relate  to  each  simulation  function. 

T'n,  ■  Tiirtjct 

The  nominal  local  ion  of  the  target  is  defined  in  terms  of  its  range  (R)  and 
altitude  above  a  spherical  earth  (HT).  Its  velocity  is  incorporated  in  the 


so 


simulation  by  defining  in  which  hopple r  filter  (KFU.T)  the  target  appears 
(KFJhT  *  1  for  dc),  The  average  radar  cross  section  (RCS)  of  the  target  is 
also  specified,  as  well  as  a  fluctuation  parameter  (KSVIfcR).  The  fluctuation 
is  assumed  to  be  slow  (pulso-to-pulse  correlated  but  scnn-to-scan  independent), 
and  the  following  cases  can  be  handled: 

KSWKR  Case 

0  nonfluctuating 

1  Swerling  Case  1 

2  Swerling  Case  3 

>0  chi-square  with  2*KSWER  degrees 

of  f  reedom 

The  target  position  as  described  above  is  nominal;  for  each  statistical  replication 
the  actual  position  is  varied  randomly  throughout  the  azimuth  beam,  the  range  gate, 
and  the  Doppler  filter.  The  beam  scanning  and  straddling  losses  are  thus  incorpo¬ 
ral  kI  in  the  simulation. 

Sp.il  lallv  *» inhomogeneous  ground  clutter  can  be  simulated.  The  model  used  in 

thi  •  program  is  a  block  cor  i  e  L . 1 1  ed  one,  in  which  a  value  of  o  (the  bacRscatter 

o 

>  fit  it'iil)  if.  generated  for  one  square  of  size  RCOR  from  a  statistical  distri- 
but  ion  l  hat  can  he  log -normal  (SVI»0.)  or  Weibull  (SVP<().).  The  for  adjacent 
squares  in  the  format  of  Figure  12  is  generated  independently,  but  within  any 
square  the  is  constant  (homogeneous)  with  an  assumed  Rayleigh  amplitude  dis¬ 
tribution  applying  to  each  elemental  area  within  the  square.  The  global  average 
ci  is  defined  as  SICM  on  input.  For  homogeneous  ground  clutter  set  SVV^O. 

Two  models  for  the  clutter  spectrum  (not  including  the  scan  modulation)  are 
fu»  l ud^d .  They  are  a  Gaussian  shape  (IT-0)  and  a  general  shape  (IT^O)  of 

*;<!)  -  ' 

1  +  07l())" 

win  ir  n  IT  (the  main  plug  taut  restricts  the  choice  to  IT  =  ’.!  or  3).  The  half- 
powet  spectral  width  relative  to  the  J*KF  is  defined  as  KWTR.  The  above  descrip¬ 
tion  lor  the  flutter  spectrum  is  designated  as  the  ”ac  component.”  There  is  also 
a  M'b-  c  «•;{«(  »o»i*  'in  the  powet  of  which  is  IX!  AG.  times  the  power  in  the  ac  component. 
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The  user  must  determine  if  there  will  be  any  ground  clutter  at  all  because 
of  terrain  masking,  and  whether  there  will  be  any  enhancement  of  short  range 
clutter  in  range- ambiguous  situations  (R‘»C/(2.*PRF) *  where  C  is  the  propagation 
velocity) . 

Rai >i /Chaf f  Chit t ev 

The  volumetric  clutter  is  assumed  to  be  homogeneously  distributed  in  a  layer 
between  two  altitudes  (Hi  and  112,  112  >  HI)  relative  to  a  spherical  earth.  The  re¬ 
flectivity,  or  backscatter  cross  section  per  unit  volume,  is  defined  as  REF.  For 
rain  clutter  set  H1“0.  A  Gaussian  shape  is  assumed  for  the  fluctuation  spectrim. 
Wind 

The  spectral  width  of  all  types  of  clutter,  and  the  mean  Doppler  of  volu¬ 
metric  clutter,  are  functions  of  the  wind  speed  (VW)  and  direction  with  respect 
to  the  auteima  beam  axis  (ANCl.). 

The  dyrtrrn 

The  radar  system  is  defined  by  the  following  parameters: 

PT  =  peak  transmit  power 

C0  *  peak  antenna  gain,  one-way 

Wl.  -  wavelength 

BW  pulse  (resolution)  bandwidth 
PRF  -  pulse  repetition  frequency 
IT.  ~  product  of  noise  figure  and  rf  losses  (L>0) 

AX1DB  -  half-power  width  of  azimuth  beam,  one  way 
l TYPE  -  type  of  scan  (-1  for  step,  =2  for  continuous) 

Al'.UWKI.  *  azimuth  angle  through  which  the  antenna  steps  (1TYPE=1)  or  scans 
( IT Y PL- 2)  from  dwell  to  dwell* 

PSAT  -■  receiver  saturation  level,  after  pulse  cancellation  but  before 
Doppler  filtering 

NGOil  i:  number  of  pulses  cohctcnLly  integrated 
NNGOH  *  number  of  pulse  groups  noncoherent ly  integrated  at  same  frequency 
NGNC1  -  number  of  stages  of  pulse  cancellation 

NR  =  iHDiiber  of  lange  gates  simulated  (=1  for  conventional  threshold 
detection,  >1  but  odd  {or  CFAR  processing — the  number  of  CFAR 
ref ei erne  cells  is  NR  1) 

I -AW  "  envelope  detection  law  ("l  for  linear,  -2  for  square  law) 

A 

a  dwell  eons  I  slo  «•!  Nmil*NNCOlI  INCNCI.  pul.es. 
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The  processing  stages  in  the  receiver  are  assumed  to  be  in  the  lol lowing  order, 
beginning  at  the  front  end: 

•  range  gating 

•  pulse  cancellation 

•  saturation 

•  Doppler  filtering 

•  envelope  detection 

•  threshold  detection  or  CFAR  processing 

The  resultant  clutter  spectrum  for  the  scanning  antenna  (ITYPE*2)  is  modeled  on 
the  basis  of  the  Gaussian  shape  for  the  beam  and  input  spectrum. 

Simulation 

The  number  of  statistical  replications  under  the  same  nominal  conditions  is 
specified  as  NREP.  For  false  alarm  analysis  set  RCS^O.,  and  importance  sampling 
can  be  invoked  with  KSW=*l  (conventional  sampling  will  be  implemented  with  KSW*0). 
The  distortion  parameter  for  importance  s.impling  is  XM,  which  should  be  established 
by  trial  and  error  (loo  small  of  a  value  will  result  in  infrequent  false  alarms, 
while  too  large  of  a  value  will  cause  underflow/overflow  or  other  diagnostics; 

XM=2,  Is  a  reasonable  first  guess).  As  set  by  a  DATA  statement  the  number  of 
azimuth  samples  simulated  for  clutter  is  NA“21,  and  the  spacing  of  the  samples 
is  computed  in  MAIN  as  DA*.  2*AV,3DB.  The  user  can  decrease  the  running  time  by 
choosing  NA»ll  and  DA=.4*AZ3DB,  or  NA*9  and  DA=.5*A23DB,  to  encompass  +  2  beam- 
widths  of  the  azimuth  mainbeam.  Further  out  sidelobes  can  be  accommodated  by 
choosing  NA  and  DA  accordingly. 

Debug  printouts  can  bo  obtained  by  setting  IDBC-O.  In  addition,  it  is  pos¬ 
sible  to  obtain  a  printout  of  the  signal  spectrum  in  the  first  range  gate  (or  any 
range  gate  with  suitable  program  modifications). 

2  1  Nn-RPRETI_N/;  HIE  RESULTS 

For  conventional  threshold  detection  (non-CFAR)  the  output  of  the  Doppler 
filter(s)  under  lest  is  normalized  by  the  Input  noise  power  divided  by  NCOH, 
which  Is  the  approximate  (receiver)  noise  power  in  a  Doppler  filter  on  output 
(It  would  be  the  noise  power  If  the  filter  weights  were  uniform).  These  detected 
outputs  are  then  accumulated  in  a  histogram,  and  the  histogram  is  integrated  to 
obtain  the  cumulative  distribution  function  versus  relative  threshold  setting. 

For  CFAR  processing  the  output  of  the  Doppler  filter  under  test  is  normalized  by 
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ihi'  average  of  the  NK-1  range  gated  outputs  Conning  the  CFAR  reference  (half  of 
tile  gales  are  in  front  of  the  range  gate  under  test,  and  half  are  behind  It). 

In  order  to  determine  what  the  threshold  setting  should  be  set  RCS=0.  to 
get  the  Simulated  false  alarm  probability  versus  threshold  setting.  Usually,  the 
false  alarm  probability  will  be  given.  Then  with  RCS>0.  we  can  read  out  the 
probability  of  detection  that  corresponds  Lo  the  desired  threshold  setting. 

3.3  EXAMPLES 

If  NNCOII-1,  KSWER=1 ,  LAW =2 ,  and  NR=1,  all  signals  appearing  in  a  non-clutter 
region  of  the  Doppler  spectrum  will  be  Rayleigh  amplitude  and  the  probability  of 
exceeding  a  threshold  that  is  normalized  by  the  average  output  noise  power  is 
given  by 


P(T)  = 


-T/(l  +  SNR) 


(114) 


where  T  Is  the  normalized  threshold  and  SNR  is  the  output  signal-to-noise  ratio. 

For  the  first  example  we  will  examine  the  false  alarm  performance  with  im¬ 
portance  sampling.  The  simulated  parameter  values  that  deviate  from  the  built- 
in  data  values  are  RCS=0. ,  8 100=0. ,  NC0H=1 ,  KKILT=1 ,  NREP=1000,  KSW-1 ,  and 
XM=13.8.  The  simulation  output  is  given  in  Table  2.  Since  SNR*0  in  (114)  we 
can  easily  compute  what  the  theoretical  results  should  be.  For  example, 

P(8.0)  =  .000315,  P(10.0)  =  .0000454,  and  P(12.0)  -  .00000614;  all  values  fall 
within  10%  of  the  corresponding  values  in  Table  2. 

For  the  second  example  we  will  repeal  the  above  conditions  except  we  will 
utilize  foul  Doppler  litters  (NCOII-4,  KFII.T=3).  The  simulation  output  is  given 
in  Table  3.  Because  ot  Hamming  Weighting  in  the  Doppler  filtering  process,  the 
noise  bandwidth  of  a  Doppler  filter  is  about  25%  larger  than  the  PRF  divided  by 
tire  uumbirr  of  fillers  (N(.i)J)).  if  we  thus  multiply  the  values  of  T  in  Table  3 
by  0.80  prior  lo  the  use  of  (114),  we  will  got  good  statistical  agreement. 

For  the  third  example,  we  add  a  target  (RCS*=4.)  so  that  the  input  signal- 
t»i  noise  ratio  is  a  faelot  NCOH  higher  (4  *  8.0  =  32.0),  reduced  by  the  following 
efforts 

Antenna  scan  .69  (-1.6  dB) 

Dopplci  filler  straddle  .90  (-0.5  dB) 

Range  gate  straddle  .58  (-2.3  dB)_ 

Total  .36  (-4.4  dB) 
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I  A I  I  U  iilluvi  ,  r  1  ,t.  >•)/.  ,V|.  )-■  ,SV1'  ,  |  I  .  CL  At  .  I  » 1 1  v.r'  ,  0.  I 
1.  At.I  ASUV  I  ST'  T  ,  «  •  ,Yk',  ni,  I 
'.'ll  A  >  Ilf  1  •  ?(.'  J  ,  •  •  Y  1  .  Y  |  I 

1  >  l  1 1  :  i-  .1  •  .0 1  :  Tl.  S'. 

C.M  I  f-  -1  1  I  Y<,M  I  .M-  Y-- .'Ill  1,  ItM.I  tl  CL  HI  I*’ -(,»  I 

l  At  l  i  Ml  I  Y  :  ,  M  t,  1  ...  «.l<  V  |  , /.IHM  71*  U  I  |l  ri  Cl  f.’-Cf.'  1 

c.  him  i..)s  _ .  i :  ....  i  ',c  i .  i,. 

I  'tl  lllill-'  |Y',YI.*..‘tl,.l,t  I  .  i.  ,t.  •  .  1  h  !  v  Vl  .!  n  ]  -  ) 
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i  M  L 

l.'illMMill  ,  >J  | 

i  Ail 

Ai'j't  Mil  .  «  !.»  1  .  M  1 

1 1  l  1 1  ;•(, .  L  ■  .  ’  1  v.  I',  i<. 

L  ••  l  1 

H  '  V  1  1  W  »  *1  ')  *  %  —  S  «  tj i « 

>y  ,  r  OliliU \  YU  \  (J 

U  11  1  IK  -  »  N  1 

l  -U.I 

»**  '*  1  1  ^  I  *  K  1  U  1  «-<m(i^ 

y  ; ,  /'ifa  i'mjt  1 1 

UUHr'-^N  » 

Jr  1  >  < 

IMr  .('.il  CO  111  1H 

W  l 

.‘■'.I/UIIM-'IM  IC.I-l. 

1 

L  M  l 

lUfl  lIVK.rl.'Ii'Crt.KSKIW.i'l  lLl.A/i.ljA/l 

<  M  l 

«MIMMfll,»H,U,»Kl 

f  AL  l 

11  l  ll  l-G.ll  .Cl  f.J  Kl  Ii* 

LAlI 

p.:miij,m 

YMvOM'UPUl  Lf 

I  AX  Cl  T  1 

CALL 

IV  Ml  (  Y  I  ,  m  T  ,  ,.l 

V 1 ,/OfCL TPUl  Cl 

YAkCfT  1 

t  *t  \  i 

PI- A  s-*rr.  UAL  AM 

1  YSIS  C  N  IK  Mill  P/..SGE  Pit-C 

■  M  tAL; 

>.v  1  1  1- V>  ,'U  ,  Y.:  i 

L  U  L 

>x  1  1  l-NS  .11.  Ml  I 

LAL  t 

y  “  j  J  imv.  s .  >p,  yk  i 

l  M  i 

'“111  APCLS  .  X  I  .  Y  1  1 

1  Ai  l 

SPCI-YIYB.Y 1 .MULSMS 

.  .OH  .  Kf.  T  1 

I..M  l 

AMCI  AS,  YB  .  S  «S  1 

111  At 

PCMCl  SS  TlO  j  1 C N A L 

CA.  I  P<M)(  (  M  «  a  <  V  t  <  *  UK  <UL  A*  .M.CIJn,  ACNCL,  PSA  T  1 
NX-  Kc'-x  |  »  1 

A  "  l  • 

IltKSfc.ti.n  Vi  -  >  •  NMi»  )  <|  XX  (  -  <  1  I.  />  "  >  ) 

P  Oi.l  T  v.  :  ’  r 

1)‘  (  KMi'k 

<  a 1  i  i’sii'MU-'t  -  >  / l  .vi  i-  n ,  w l 

AG  i  i  \  T  1  MJ * 

*  0  (.  U.'  1  I’D! 

i  i'r.  i  \t  ii.i  :i  i  |*",1  xam.i  oikO 

. *  .  -v-  .  .  /  /  \ .  M 

l.  .‘.I  I.  Cvir  <-NS  .  SMWV.Y  .  ■>  » 

l  .'il  M.,l(  I.\ui)iin  01  IM  Ull  I 

(  -’>1M  1,1.1  IHi'UHUll'  MAtlSIJL) 

C  -u  L  htlliMHI 
MM  110 
ii :  >■  i.  J '  1  •  M 
i  •  1 1  •  t 

p  w:.u » l  .-r  1 1  j  i 

!’lIM  1  1  1  .  t  .  h!'," 

(.(•  !  i .» n  vj* 

tr  •> 

I  •<  MM  1  UU  \Si;>  ,  «.S,iv 
Ml  P 

M.  I  ■  -  Mill/Wi  liv;  Milt  1:  U'l  '.Slilfi  SUI  IS1  ICW/ 

!  t  >  Mi.il  1 1 .  A  I  /  i  0*  1M1  ,  ?»'.t-OA  1  A  ,ct  .  T  /  1 

i:i  t  ;  <  “ ;  1 1 '  i  u  <• , ;  :  •. . "  i 

II.'  ,  :  *■  V  T  (  /  v  /  /  SM  '••«*  1  1 1.  M.I  U-1'Ul.i,  MJl.rlT  1>  M  C  A  Y  1  V  I  / 

1  ■  ••  '  SI  . . h./HMI 

i  I 


.1 
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S  111  INI  I  At  ;,l  I  (  X  1  .  1  -M  .  k.kf  1L  1  ,A/  1,:>AM 

C  IN  lf*S  SIMNl.li  IJM  ll-l  Ct"Vlf  <  V1I>M)  SIC’FAL  FOR  Tilt  lARGIT  IS  tULAJf 
C  till  ARIAY-PAIR  (XK,X|).  1  HI  1AHGI  I  CUPr-Lt*  IS  RANCC**  HllHlK  IM 

l  Mill  1  XUC:>l>  kllll  (  1 1  NCCH.Ct.il  (W  X  I  T  H  l  S  trtf  Piff  (IF  NCGF-ll. 

C  tit  A/l*!0lH  HI  A*  SaFFPS  HAST  fHf  ItHI  1  UUk 1 NG  lHE  U»E  L  L  LF  NPLLS 
C  PCLSFS  HtGINMNC  At  A/1  »11H  AN  A/ I  “UTH '  1  NCR  E  HE  Nt  OF  CAi  FOR  EACH 

C  PCISI  .  1  HI  I  CS  PAKA*-(  I£HS  ARE . 

C 

C  TPPx  ■  lAKClt  AVlRAGt  X I  C  f  1  VI D  P|)wtR  (  •  R  SCL  «k  C  S«l  LC  A  I N*  *2  » 

C  KS«IK  •  SnfkLlNG  CHI  -SCLAV  I  PAl'.\l*tH°  (SfE  REF  .31 

C 

C  fl-t  F  (  k  r  A  t  (If  IHf  XF.XF-ARRAYS  IS  AS  IF  THEY  xMl  C  J  **E  N  S 1  CN  f  C 
C 

C  XxINPULS.NR)  X  I  INPULS. NRI 

c 

C  IXI  AS  I  .1*.  >  C  I  1  I  ,  x  I  I  1  I 
C  «“|N  /  t  I)  “  l  /  M  .Mi't/i  ..NA.NCIIH 
c;'“*'CN  /i  hoc/  I  ::H(, 

C  A  t  A  I  it  OH  I  /ri./sJI  -IS  1/ 

PIU  =  1PIIH 

I  It  *  Ski  ■!  .1  -.0  I  0"  If!  // 
pxOJ* I • 

Cl.  /O  J*  .  .  X  Sk  I  r 

p«i:t:*PR (’.>»« am  to.  > 

,'U  C  t  I  Nut 

PI  »  ’  1  Pli.  »  I  -  AC  (I,.  I  PSl;i)|  | 

?.'  A-scRiin.ii 
(•  -  .  S*RANI  I  i!  .  ) 

t  (A  I  It  f  *  *■'  A  N  f  Ki  .  )-|  .  S  l/f  CIA  1  I  NCDH  ) 

ill  l.R-  1  I  //  )  St  1 

If  I  li'l'o.t  l  .01  P.|M  |OO.A,l-,l 

(.At  L  >1  t  r  C  -  Nk  *M’l!l  SfC.tXkl 

I  AH  >'*1'  t  -N>  *NI*Cl  S  .<!.  .  X  I  I 

A  /  A  /  1 

i  i,  ,H  K  “  1  .Sf'Ul  S 
ARC  -  I  FOP  I  *K  *1 
S R-f PS l AX  Cl 
S  I  S  I  M  A  0  C  I 
ARi>  -  A*A/..A  1M  Af  | 
x  S'  I  1  I  >  A*P»  I  !.  -!•  I*  SC 
<1(1  I  »  A  •*••>♦  I  I  .  -HI  *  SI 
I  >  1 1. is . l !  .  l  I  c;  id  , 

XP  I  l  .NPC,  s  I  tA  «P»|-»SK 
X  I  1  l  tNPCI  S  I  *A.“i>Y|-»S  | 

/*>  l  -I  t  1 

r.f  a/.;ja/ 

111  Cl  N  1  I  l  1! 

-Ml:  if- 

1011  Ft  »**Alf ////1H  IA-l.,1  A.F-.I . I  l/.At/t  l/.FI 

I  M 


.1 


;o 


SUBROUTINE  CLUTTRIXR.XI.C^CH.PCOfi, SVT, IT,OCAC,FHTR,POTI'» 

r 

C  IN  TH5  SUBROUTINE  WE  GENERATE  THE  CCMFLEX  VIOEO  SIGNAL  FC=  CllTTEF 
C.  WITH  A  STfF  SCANNING  ANTENNA.  THE  CLUTTER  SIGNAL  IS  CLTPUT  TC  THE 
C  ARRAY  FAIR  (XR.XII.  THE  PARAMETERS  ARE..... 

C 

C  CPOH  =  CLUTTER  RCS  PER  UNIT  AREA  ON  THE  GROUND  SCALEO  TC  RECEIVEC 
C  POWER  «-SIG0*FSCL1 

0  R  *  RANGE  Of  INTEREST 

f  OR  -  RANGE  RESCLUTICN  CELL  SI7£ 

r  OA  -  A7IHUTH  SAMPLE  SPACING 

C  NA  -  NUMBER  CF  AZIMUTH  SAMPLES 

C  RCCR  -  CLUTTER  CORRELATION  DISTANCE 

C  SWP  =  SFATIAL  VARIABILITY  PARAMETER 

r  roTR  r  mean  dofpler  of  clutter  relative  to  the  fff 

C  PHtr  -  oOrFLER  WICTH  OF  CLUTTER  RELATIVE  TC  THE  FRF 

C  IT  i  PARAMETER  SPECIFYING  SPECTRAL  SHAPE  «SEE  SLR .  FANSECt 

C  DCAC  -  RATIO  OF  OC  POWER  TC  AC  F CWE R 

C 

r  IN  ADDITION  n  is  the  number  of  spectral  samples  used  in  the  GENERATION 
C  OF  THE  r.ORRfLATEO  RANCCM  SEQUENCES.  IT  SHCLLT  ME  AT  LEAST  AS  LARGE 
C  AS  M  ULSU.  7FWTR. 

C 

C  THERE  ARr  TWO  GENERAL  CASES  OF  INTEREST . 

0 

C  GRCLNO  CLUTTER  RAIN/CHAFF  CLLTTE' 

0  - - -  - ........... - .... 

C.  SVP.NE.C.  INGNHOMGCE  NE  OUST  SVF-C.  tHOMCGENEOLU 

f  FCTP-0.  1CTRi|Z»VWFWL» »ccsiangu 

0  nCAG.CT.O.  DCTC-0. 

r  IT  -  C 

0 

C  IF  THE  EIFVATICN  REAM  COES  NOT  HAVE  PCAK  GAIN  (N  THE  GROUND  TT  THE 
C  RANGE  CF  IMPOST  THEN  CPOW  MUST  INCLUCE  THE  ELSVATICN  REAM  WEIGHTING. 

C 

f  THE  FORMAT  OF  THE  XR ,  XI- AP. SA Y S  IS  AS  IF  THEY  WERE  Cl N  ENSI CNE  t 
r 

C  XRENPLtS.NM  XT  (NFLLS.NPT 

0 

nuFNSiOK  xMii.xun 

COMMON  /r.OMl/  NR,NPULS,NA,NCOH,R,0K,CA 
COMMON  /ISAM/  KSW ,ISW 
COMMON  / OR U G 7  lORG 
OATA  T WOP  T/6.7AT1FS3/ 

NO- (  NP  -*1 1  /E 
F-rPON*P*OR»CA 
HRrSANr  »0.» 

MR“  - 1 
11  =  1 

GO  '.0  TM.NR 
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IMMfi.EC.IFIXIRRt  I  GO  10  25 
SUH=0 • 

A7  =  «.5MNA-11»0A 
AA^RANT (3  .  | 

MA~0 

mSPTVARlSVF) 
no  20  J~1  .NA 

IFIHA.EC.  1PIXIAAH  GO  TO  l*i 
W=SPTVAR«SVP| 

R.A=AA 

1=  SUM=SUR+P*N*A7GATN ( A7> * *2 
A7=A7»TA 
AA=AA+R*CA/KCGF 
20  CONTINUE 
MR-RF 

25  PAV=SUH 

IFIIO°G.GT.O »  FR I N  T  IOC, I, PAW 
ISW  =  0 

IFIKSM.EO.l  .ANO.I.EO.NCI  ISN  =  1 

CALI  RAKSECtXMLll.XIUtl  ,NPULS,  FM 1  R ,  It  ,OC  AC  ,F  AW  ) 
Ir<rOTR.eO.O.I  GO  TO  35 
L  -L 1 

no  30  *=1.NPULS 
ARGrTNCFI*FOTR»lK-l» 

C=COS IARG ) 

S=SIMAPG» 
fR  =  C*XRa»-S*XT»L1 
XI-C*XTIL» ♦POXFTLI 
X  R  (L  I  =  VR 
XHL)  =  YI 
1=1*1 

30  CONTINUE 
75  RR=RR*rfc/RCOR 
L  1  r  L 1 +NFUL  5 
40  CONTINUE 
RETURN 

100  FORMAT (///2 IN  CLUTTR  I.FAV . I5.E12.4T 

E  NO 


r>r>n^^>r>nr>nn 
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SUFROMTINE  NOISE<Xfi,Xl,PN> 

IN  THE  StJFRCUllNE  THE  THERMAL  NOISE  I!  generated  in  tne  aerty-paif 
<X»,XIt. 

PN  -  AVERAGE  NOISE  POWER 

THE  FORMAT  Cr  The  Xfi.Xl-ARFAYS  IS  AS  IF  THEY  HERE  CI*ENSIGNEC 
XRINPLLS.NR1  XII NFLLS  *NFY 

i  ihenficn  xRin.xnn 
COMMON  /ISAM/  KSW.ISH 
FCFMON  /r.CKl/  NR.NFULS 
Nf.=  CNR  *11  /2 
L-1 

TO  ?0  1=1 , NR 

1  sw  =  o 

IFIKSH.EO.l.AND.I.tD.NO  ISW=1 
'  0  70  K=1 ,NPULS 
I  ALL  GAUSS!  (X®U  I  .XT  It  I  ,PK» 

1  -  L  *  1 

?C  ONTINUE 
ETllcN 
'  NO 


f 


SUBROUTINE  FfiOCfS<XF,XI,Kl,K2,LAN,NNCCK,NCNCL,FSAT» 


r 

r  IN  TV  IS  SUPKGUTINE  THE  SIGNAL  FRCCESSIRG  IK  THE  RECEIVER  IS  1HFLE- 
C  M.LNTrO  ON  THE  ARRAY-FAIR  IXR.XTT.  THt  ARGUMENTS  ARE.,... 

r 

C  K1  -  FIRST  FILTER  FOR  OCPPlfP  SEARCH 

r  K?  -  LAST  FILTER  FOR  nCRPLER  SEARCH 

C  LAH  =  LAM  OF  FIPST  OETECTOR 

C  NNCOH  -  NUMBER  CF  FLLSE  GPCL°S  NONCCHERFNTLY  INTEGRATED 

C  NCNCL  =  STAGES  OF  CANCELLATION 

C  FRAT  =  SATURATION  FOKER  LEVEL 

r 

f  THE  FOLLCHTNr,  CASES  CAN  BE  HANDLED . 

r 

C  FUISE  CANCFLLATION  HTI  NITH  BtNARX  HEIGHTS 
C 

C  NCNCL  -  STAGES  OF  CANCELLATION  tNCNCLM  =  NUMBER  OF  FLLSEST 

F  NCNCL  =  0  CCRRESFGNOS  TO  NO  PLLSE  CANCELLATION 

C 

C  FORflEF  FILTERING  via  F r t 

c 

C  NCCH  =  SI7E  OF  FFT 

C  NCOF  r  1  COR CF SFGN CS  TO  NO  OOPFLER  fxltFFING 

r 

C  SATURATION  P  CTwt  £ N  FULSE  C  ANCEl L  A II C  N  AND  DOPPLER  FILTF:ING 

r 

r  FSAl  -  SATLRATION  fower  level 

r  FSAT  =  1.EB9  CCRRFSPONtS  TO  LINEAR  PROCESSING  INO  SATLRAT ICM 

C 

C  LAH  OF  FIRST  CETfCTOF 

C 

r  LAN  -  1  CGRRESFCNCS  TO  LINEAR  DETECTOR 

C  LAH  =  2  CORRESPONDS  TC  SOUARE-LAH  DETEETOR 

f 

C  NONCOHERENT  INTEGRATION  AT  SAME  CENTER  FREOLENCT 
f 

r  NNCCH  =  NLHF.FR  CF  FLLSE  GROUPS  N ONCCKER ENTL V  IN'TEGF  AT-f  C 

r  NNCOH  I  1  CORRESPONDS  TO  NO  NCNO.OHERFNT  INTEGRATION 

C 

r  f  F  Afi  THM  SHCLO  RECLLAT  ICN 

C 

C  NH  -  TOTAL  NUMBER  OF  RANGE  CELLS  FROCESSEC  INCILCING  TARGET 

C  NR-1  -  NUMBER  CF  GEAR  REFERENCE  CELLS 

r  NF  -  1  CORRESFCNCS  TO  CONVENTIONAL  THRESHGLC  PROCESSING 

f 

r  THF  INPUT  complex  SIGNAL  IS  in  THE  array-pair  TXF.XII,  THE  format  CE 
f  WHICH  is  AS  IF  THE  DIMENSION  HERE  XR I N PULS ,NKT ,  X  I IKPULS ,NR » . 

C 

r  the  noncohff.fmly  integrated  output  cf  filters  mi  through  K2  affears 

r  IN  .'FRAY  XF  IN  SAHFLES  1  THROUGH  K2-K1»1.  IF  NG  DOFFLER  FILTERING 


?A 


r  VTA  F  r  T  rs  IHPLENCNTEC  THEN  NCGH=H  .  USUALL*  I  f  A  I  ARC  Et  IS 

C  ‘RESENT  THEN  Kl^x?  Hitt  CORRESPOND  TC  THE  COPPLEF  SAHFLE  IK  HHICH  THE 
r.  target  has  placed,  if  no  target  is  prfsent  then  there  ARE  SFTEPAt 
C  POSSIBILITIES..... 
r 

C  1.  NO  clutter 

r 

r  SET  HI -1 .  K?=NCOH 

C 

r  2.  r.RCUKC  CLUTTER  CKLV  WITH  HAHHIKG  FILTER  HEIGHTS 
r 

C  SET  K1-.I,  K2=NCOH-l 

C 

C  T.  GPOUKC  pLUS  HEATHER  CLUTTEP 

C 

C  SET  K1  -  FIRST  DOFPL E°  SAHFLE  NGT  FLAKFET 

C  <2  -  LAST  C)OPnLFR  SAHPLE  NOT  P|  ANKFO 

C 

C  for  CFAR  PROCESSING  THE  RANGE  SAHPLE  UNDER  TFS  T  IS  ALWAYS  THE  TENTER  , 

C  SAHPLE. 

C 

DIMENSION  xRm.XItll 

riMENSioN  ;tiebi 

CGH-ON  /r on/  NH.NFULS.NA.NCCH 

COMMON  /oouc,/  innr. 

IF  (NFULS.NE  .  KNCGM  »KCOH  .NCNCL  T  STC®  1 

IF  (MCGTNP.E* .NE.1T  ST CF  Z 

!F(NrOK.£C,l. ANO.Kl.NC.il  STOP  3 

IFINCOH.EO. l.ANC.KE.NE .1)  STOP  A 

SREFrKC.j 

MP-NPUl S 
'  C  SO  I  - 1 ♦ K R 
1-  f  1-1  I  •NPULS.l 

IFCIDBG.Lt.a.OR.I.GT.ICRC.T  F.G  TG  2? 

CALL  FRNT<XR<L1I .NPULS.-N.GH  XR,20HINRUT  TC  FRCCES  I 

CALL  FRNTTXl  (111 , NFULS.-A.fiH  XI «  2QHTNPUT  TO  PRCCES  » 

22  I  F  (  NCNTL*  E  C  .0  I  GO  TO  2N 
C  MERC  HE  INCLEMENT  PULSE  CANCELATION 

TO  2S  H  =  I,  NC  NCI  « 

l-Ll 

'.F.-NP-i 

00  2  0  Hr  1  ,KF 

XR(LI-XR(ll-XfiU*l» 

XITLTrXIlt )-XI 1L‘ 1 ) 

1-1*1 

20  FENTINOF 
2S  CONTINUE 

IFITDPC.LE.O.OR. T.GT.IDRGI  GC  TO  2 f 

CALL  FHNTIXRfLlt .NF.-A.6H  XR , 20 H A F TER  CULSF  CANCEL  1 

FALL  Ft  NTTXI  (LIT  ,NF, -A, F.H  XI  ,  20  H  A  F  TE  R  PULSE  CANCEL  »  j 


?{■  IF  (PSAT.GT.l.ERiil  GO  TG  3? 

C  HERE  WE  IHRtEHFNT  SATURATION 
L=L1 

HO  30  *  =  1 » NP 
I=TRlLl**2«XItLI **2 
IF(P.LE.PSATI  GO  TO  30 
A=SOKT (PSA  T /FT 
XR(LI=A*XR (LI 
<1 (L ) -A»XT <L» 

L  =L  ♦  1 

*0  rONTINUt 

IF (ICPG.L F. O.OR. I.GT.IGRGI  GG  TO  33 

CALL  FRNTCXR  (Lll  ,  NF.-4  ,6H  XE. ,  2DFAFTER  SATURATION  I 

FALL  FRNdXKLlI  ,NP,-4,6H  X 1 , 20 F A FTER  S ATUR AT ICN  t 

r  HCR E  HE  IHFLEMENT  COHERENT  FILTERING  AND  NONCOHERENT  INTEGRATION 
33  L  =1 1 

CALL  XFIT(-NC0H,D.,7) 

CG  40  K-l « KNCOH 

CALL  FfHFLT (XR1LI . XI 1LI ,L AH, NCNCLI 
CALL  ATUF (NCOH.XR (LI ,7,71 
L-L+NCCH 
40  “ONTINUE 

IF  (  ICFC..LE.  O.OF.  I  .CT.IOF.'.!  GO  TO  42 

(  ALL  FRN7  (  XR  (L  II  ,NF,-4,6H  XP , 20 F A F TER  C CH  F ILTER  I 

42  .  ALL  XMT  (NCGH,  7,  XRIL1 1 1 

SO  i'CNTIME 

<KrK?-Kl*l 

L1=F1 

11  =  1 

~C  55  1=1, NR 

ALL  XFTT (KK.XFl Lll ,XR (Hill 
\.1  =  L1fNFULS 
11=M1 ♦  KK 
55  '.ONTINUE 

IFdORG.Lr.0l  GO  TO  57 

ALL  Fcst (Xfi,NF»KK,-4,6H  XR,20F6FTEfi  OF  Trn ICN  I 

57  IF (NR.FC, II  GO  TO  75 

r  here  hi  ipflekcnt  ctar  prccesstng 

•■RG  -  (NF.l  l  /? 

C  70  Krl.KK 
*  F 

SUF=C. 

C  F> C  1  =  1, NR 

IFII.NL.NPCI  Sl“.=  SLH  4XR  a  I 
i.=L  *KK 
50  'ONTINUE 

.  -  INRC-ll *KK«K 
*0 ( Kl =NR£F*  XR (LI /SUM 
70  'ONTINUE 

IFdDPG.Lf, 01  GO  TC  75 

CALL  FF NT ( XR  ,KK,4 ,6H  XR,20HAFTCR  CFAR  » 

75  4  F  TURN 
r  NO 
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SljPROUTINF  OGHFLUXR.XI  .LAW.NCNCH 

r 

r  IN  THIS  SLFROUTINE  WE  IMPLEMENT  COHERENT  FILTERING  CR  THE  CCNFLE  > -Vlt  EC 
C  SIGNAL  IN  THE  ARRAY-PAIR  ( XR , X I )  OF  LENGTH  NCOH.  THE  OETECTEC  OUTPUT 
C  APPEARS  IN  ARRAY  XR  OF  LENGTH  NCOH. 

C 

C  NCOH  =  SI7C  OF  FFT  IN  COFFLER  FILTER  EANK 

C  NCNCL  =  NUNEEF  GF  STAGES  OF  PULSE  CANCELLATION 

C  LAW  -  ENVELOPE  OETEl'.TION  LAV  1 1  =  L 1NE  AP ,  ? = SOU AK E -L AW> 

C 

C  PUL  St  WEIGHTING  IS  APPLIED  PRICR  TO  CCFPLER  FILTERING.  THE  RFIGHTING 
C  FUNCTION  IS  COSINE  ON  A  FECESTAL,  WHERE 
C 

C  AlPHA  -  1.  FCR  UNIFORM  WEIGHTING 
C  =  .OP  PGR  HAMHING  WFIGHTINF. 

C 

1  tmension  mmi.xmi 

(•  IHENS TON  C  *12 fl»  ,W  (12fl) 
f.  OH  MON  /COM1/  NR.NFULS.NA ,NCCH 
■ATA  ALPHA/. 08/ 

AT  A  f  (11/1./ 

'  ATA  M,N?/0»0/.nI/3.1<*lGPEH5/ 

•fincoh.eo.it  GO  TO  IT 
F INI .£C.NCCH.ANP.N?.EC . NCNCL T  GC  TO  IE 
M ■ NCOH 
NCNCL 

ALL  WEIGHT  CR.NCOR.ALFRAT 
(1  )  -  1  . 

0  1?  Kr’.NCCH 

"  -=STN|P1*  f  K  - 1  T  /  Ft  CAT  TNCOHM**  (3fNCNCLT 

13  •-ONTINIJF 

IE  ALL  Ft-fO  TNCCH.W.XR.XP) 

ALL  PRG!)tNCOH,W,XI,XI) 

AIL  FFT2IXR.XI.NCCH, -IT 
17  ->0  ? S  K  =  1,NCCH 

(R (KT  =YR(K> ••2*XI IK»  »*3 
<Rixi=yRix»/riKi 
(F (LAH.E0.1T  XR1XT -SCPT1XRTXM 
7r-  '■  F'N  T  INliF 
tTIIPN 
■Nr 


o  o  n  n  o  o  o 


FUNCTION  A7GAJNTA7F 


COMPUTES  ONE-MAY  POWER  GAIK  AT  A 7 IMtITH  ANGLE  AZ  (RAC).  GAIN  IS  NOFN 
ALIZEO  TO  A  PEAK  VALUE  OF  UNITY. 

IN  THIS  EXAMPLE  A  GAUSSIAN  IE  AM  SHAPE  IS  ASSUMED  WHERE 

A7.TPP  =  ONE-MAY  HALF-POKER  MIGTH  (RADI 

DIMENSION  G<47> 

COMMON  /COM1/  CUMMY17».AZ3r.P 
TATA  NN/07 

1FTNN.GT.C*  GO  TO  25 
NN=1 

CO  20  1=1, A2 

C  tTI=EXFI-(l.P.651M.65*II-l»  >»**21 
20  CONTINUE 

?=  A=?0.»AFSI A7I/A7SGE 
IF(A.GT.40.»  A=-40. 

1  Ar  A 
A  -  A  —  I A 
T  A- I  A* 1 

A7rATN=(l.-Al‘CUAJ*A*G(IA4l» 

F FT  URN 


r>  r>  n  c*.  o  o 
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FUNCTION  ELGAIMELJ 

COMPUTES  ONE*  WAV  POWER  CAIN  AT  ELEVATT  (N  ANCLE  El  IFACI.  CAIN  IS  NCSF 
ALI7E0  TO  A  PEAK  VALUE  OF  UNITY. 

IN  THIS  FXAMFLE  A  COSECANT-SOUARE  BEAWSMAFF  IS  ASSUFFC. 

rlMENSION  C  V  4.2 1 
f  ATA  BT2/1.5707B632/ 
fATA  ELfl/.l/ 

(  ATA  NK/G / 

1  F  (NN.GT.OI  C.C  TO  2‘ 
r  n=i 

( 0  ?0  1=1.41 
f  =  FI2MI-l)/40. 

(  III =1. 

]F(E.GT.EL0>  G(I>=  CSINIEL0»/SIMEM**2 
?0  rONTTNUE 

G142)=Gt41l 
?S  E=40.*EL/PI2 

IF1E.GT.40.Y  t  =  4  0  • 

IMF. IT. 0.)  E  =  C. 

Ir  E 

E-E-It 

IE=TE»1 

rir.AiKrn.-E>»c«iF>*E»c,nE«n 

KETIU5N 

PNC 


function  SFTVAMSVFI 

r 

C  THE  RANDOM  NUH.FERS  CORRESPONDING  TO  THE  TERRAIN  SPATIAL  V AR1 AC  IL IYY 

C  APE  GENERATED  IN  THIS  SUBROUTINE.  THERE  ARE  THREE  CASES . 

C 

C  StfP.LF.C.  HE  I  BULL  eiSTRimmCN,  -S  *»*NfI6CLL  PARAMETER 

r  SVR.EO.O.  HOMOGENEOUS  TERRAIN 

C  SVR.GT.O.  LOG-NORMAL  OISTR 1 4UTICN,  SVP=STC  OEV  CF 

C  LOG-WARIATE  CNEPERSI 

r 

C  IN  ALL  CASES  THE  MEAN  VALUE  OF  THE  RANEOH  Nt'MFER  IS  UNITY. 

r 

C  SE-  REF.l.  EOS  9.1G  AND  9.11 

r 

DATA  SS/0./ 

SPTVAR=1. 

IFTSVP.EO.D.)  RETURN 
IFISVF.GT.3.1  GC  TC  20 
IF  ISS. EG. SVF )  GO  TG  10 
SS=SVC 
A --S VF 

GA-GAMHAI A+1.1 
13  E--ALOG(FAN'F(0.l> 

SPTVAi=  =  E#*A/r.A 
RETURN 

?0  G-GALSSIOUHKYT 

SFTVAR-CXPISV',,«G-.5»SVF»  » 

RETURN 
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SUBROUT INf  RANSEQ(Xfi,XI,N,FWTR,IT,CCAC,FAV» 
r 

C  IN  THIS  SUFROU I I NE  ME  GENERATE  A  CCRREUTEC  SECUENCE  CF  GAUSSIAN 
C  RANDOM  PHASORS  IN  THE  ARRAY-PAIR  IXF.XII  CF  LENGTH  N.  THERE  ARE  TWC 
C  CASES  FOR  THE  SPECTRAL  SHAPE..... 

C 

C  ITrO  GAUSSIAN  SPECTRAL  SHAFE 

C 

C  IT.GT.O  SPECTRA^  SHAPE  IS  1/ ( 1* I2F/FW7R1 *• IT » 

r 

r  in  either  case  fwtr  is  the  ratio  or  THE  2-sioec  HALF-FCHER  SPECTRAL 
C  HIOTt  ANG  THE  PRF.  THE  ABOVE  SPECTRUM  TE  DESIGNATED  AS  THE  TC-CCMFCN 
C  ENT.  IT  IS  CENTERED  AT  DC.  THERE  IS  ALSC  A  OE-COKFCNENT  ACtEC  TO  TH 
C  ABOVE,  HHEBE  OCAC  IS  THE  RATIG  OF  OC  TC  AC  POWER.  THE  AVERAGE  POWER 
C  OF  THE  OUTPUT  SAMPLES  TS  PAV. 

C 

C  THE  RANDOM  SEQUENCE  IS  GENERATED  FT  THE  FFT- METHOC  WITH  INTERPOLATION 
C  THE  PARAMETERS  CHOSEN  BELOW  WILL  FROCUCE  A  MAXIMUM!  OF  -50tP  SPURIOUS 
r  SPECTRAL  RESPONSES  WITh  AN  ACCURACY  OF  ABOUT  2-PERCENT. 

r 

DIMENSION  XRT11  ,XI  111  ,  ART2B7I  .AH2STI  .SE2STI 

data  ff/-i./,ii/-i/,nn/-i/ 

IFTFWTR.cO.FF.RNC .IT.EC.II.AND.N.EC.NM  GO  TO  TS 

PF-FWTF 

II--IT 
NN  =  N 

IF (FhTR.LE.O.I  CO  TO  bS 
H-AMIN1 (10 . «FHTP, 1.1 
NF  =  lf 

10  NFH=NF/H-a.5 

IFTN.LE.NPHI  GO  TC  12 
NF  ~  2  *NF 

r.r  to  10 

12  IF (NF.FT.pPB)  STOP  11 
NF2 - Nr /2* 1 
l IN- 3*NF*FWTR/H 
SHIM. 

r.n  ?o  kz?,nf? 

SIM  -  G  . 

IF (K.FT.LTMI  GC  TC  IS 

IF(IT.Lt.O)  SIM  =EXP(- (1.6HS1* (W-l> • N/ EFWTK* NF * I •«2I 
IF  TIT.GT. 0»  S(K» -1./ (l.»(?.* (K-l l»H/ (FWTP»NF»» •«IT> 

IS  GTNF.O-MrMM 
20  CONTINUE 

S"M-7$UM(NF ,ST 

sin-sin  ♦tCAC*rup 
SUM-SUM*  (l  .  ♦CCAO  /r  AV 
CALL  FCCPI-NF, 1. /SUK.S.S1 
IS  IP (FWTR. I  P.0.1  GO  TO  HS 
no  (.0  H-1  ,NF 
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CALL  GAUSSI  IARIK1  ,AI«K> ,SIK» » 

AO  r GNTINUf 

CALL  FrT?(Afi,AI,AF,l» 

ARCNF*il=ARIH 

AHNF411-.HM1I 

T-O. 

DO  50  K=1,N 
I  =  T 

TT=T-I 
1  =  1*1 

XMK»=ll.-Tlt»ARtI»*TT»AMT*l» 

*TIK»=ll.-TT»*At  II»*TT»Am»l> 

T  =  T*F 

50  CONTINUF 
GO  TO  70 

65  CALL  GAUSSMAR.RT.FAVI 
CALL  XMTT»-N,BR.XF1 

CALL  Xt'in-K.EI.XI)  I 

70  PFTUPN  ! 

FNC 
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SUBROUTINE  GAUSS  I IX, T . FT 

GFNERATES  F ANCON  NUHEERS  F  CR  TfE  FOLLOWING  TWO  OISTRIFUTION  FUNCTIONS 


CALL  GALSSIIX'T.F) 

CALL  EXFI1X1 


GENERATES  FAIR  OF  GAUSSIAN  RANCCM 
PHASCR  CCKFCNENTS  CF  AVERAGE  FCWEF  *> 

GENERATES  EXFCNFNTTAl  FANOOK  VARIABLE 
WITH  UNIT  AVERAGE  POWER 


THIS  SUBROUTINE  SET  CAN  ALSC  FE  USED  WITH  IMPORTANCE  SAKPL1NC  IF  KE 
SET . 

ISW  -  1  TC  ACTIVATE  IMPORTANCE  SAMPLING 
XM  =  DISTORTION  GF  MEAN  FOHFF 

FOR  THE  FIRST  CALL  WHEN  IMPORTANCE  SAHFLING  IS  IN  EFFECT  WE  FUST  SET 

WSUN  r  0. 

NSUM  i  0 

AFTER  TMF  LAST  CALL  THE  IMPORTANCE  SAMPLING  WEIGHT  IS  GIVEN  EV 


C  W  =  <XM*»NSUMT»EXF<-«1.-1./XM»*NSUK> 

C 

COMNCN  /ISAM/  KSW , ISW ,XM, NSUM, WSUK 
DATA  7/0./ 

ITYPE-1 
GO  TC  1C 
ENTRY  FXFI 

I  Type  =  ? 

10  XMM  =  1 . 

IFlTSW.EO.lt  XFMrXM 
E^rYMM*  1-ALCGCR  ANE  »  7 »  1  ) 

X  =  E 

IFlTSW.Nf.lt  GO  TC  12 

WSUH-WSUM+e 

NSUM -NSUM 4 1 

12  IFlITVFf.GF.2t  RETURN 
F -SOFT IF»E1 
IS  A~F  ANF  (  7t 
Ar A  4A-1 . 
n-CANF  (71 

0  -  6 • B—  1  , 

A? -  A • A 

R?  =  B»P 
C  -  A  7  *  F? 

IF (C.GT.l.t  GC  TG  15 
>-f * (A2-B21/C 
«E*ATi/C 
PE  TURN 
F  NT 
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FUNCTION  GAUSSCGUKHV* 

r 

f  GENERATES  GAUSSIAN  R ANGGK  NUM6ER  CP  ZERO  KEAN.  UNIT  VARIANCE. 
C 

GATA  1/0/ 

IF  1 1 •  1.1.? 

1  CALL  GAUSSI  IA.E.2.) 

GAUSS-A 

1=1 

GO  TO  ? 

?  GAUSS1 F 
T=0 

T  RETURN 
FNf 


84 


FUNCTION  APKOCIN.A.d.C) 

C 

C  THIS  iUehOUTINi  PACK  AC  t  PRCCESScS  ARRArS.  IN  £  Vt  fcr  CASE  1AESTM  IS 

C  TFt  LENGTH  OF  ALL  ASF ATS.  THERE  A* l  SEVERAL  ENTRIES . 

C 

C  CALL  XM1T1-N.A.BI  BIKMAIll  FOR  *  =  J.N 

C 

C  CALL  XHIITN.A.ei  B(KI-A(KI  FOR  K  =  1,N 

C 

C  2S-7SUH(N,AI  /E-SUHIA(Kt)  FOR  K=1,N 

C 

C  CALL  ASUM-N.A.ti.C)  C<K»=AI 1) «e«K)  FCR  K  =  l,  N 

C 

C  CALL  AEUMN.A.a.C)  C  CK  >  -  A  (  X)  ♦  MK  >  FOR  K  •  1  •  N 

C 

C  Ct  =  OOI  (N  ,A,ei  OMSUMI  A(K)»B  IK)  )  FCR  K  =  1,N 

C 

C  Call  PROOI-N.A.t.C)  C  IK  I  =  A 1 1 )  *E  (K )  FCR  K=1«N 

C 

L  CALL  PfiOOl N.A.b.CI  C (K) -A(K) *e IK >  FCR  K=1,N 

C 

L  l=EN&Y (N.AR.Al)  t =  SU1 ( AR T K> AI  IK) «»Z )  FCR  K  =  1,N 

C 

C  CALL  PCHRIN.AR.AI, FI  P (K I = AR (K I *»2*A I C Kl **c  FCR  K=1,N 

C 

C  CALL  SRAPIN.A.Bl  A  ( K I  ANL  MX)  ARE  SNAPPED  FCR  K  =  1  ,N 

C 

C  SOMl  cXAHPLES  ARt . 

C 

C  SUM/SLM  IN, A)  NORMAL  1  2E  A- ARRAY  BY  SUM 

C  CALL  FROG  l-N.l./SUM.A.A) 

C 

C  c  cIIGy  (N.AR.Al)  NORMALIZE  C AR , A  I) -A k R A Y S 

C  ANLRM- l./SOnl <E)  BY  TOTAL  ENERGY 

t  CALL  PRCJOJ-N.ANCRN.AR.ARI 

C  CALL  PRO0I-N.ANCRM.A1.AI) 

C 

C  SLMSO- Co  1 (N.A.A)  SLM-SCUARE  CF  cLtMENTS  IN  A-ARRAY 

C 

C  N  L 1 c  IFU I . 

c 

c  t  =  LNGr (N,»x. A  1 |r  OOT IN, AR.Ak)  »OOT (N,#I, All 

(. 

o  ll-LNSION  Alll.Ull.CUl 
t M  RY  PROC 
It  IM  1  0  .  1  <’  »  1  C 
10  Nr.-  - N 
A  i.  All) 

CL  11  K -  1  , N  N 
G  O'.  >  AA'i'IK) 


ji  ctMiKue 

12  lidW 
10  LC  18  K-l ,N 
C(Ki-«(K>*B(K» 

18  CCM  IKU£ 

Hlt^K 
cMRY  cNGY 
ct-C  • 

OC  20  KM,8 
U=tt*»(KI*»Z4HIO**? 
20  CCMIMJE 
AHtiCC=tt 
kc.  TURK 
£MRY  POHR 
CC  30  K= 1 ,N 
coo  =  At*)  *•  2+eno  »*a 
30  CCM  1MJ£ 
iu  ICR* 
dlM  ZSUf- 
SUtoO* 

UC  35  K  = 1 ,N 
iCrt=SCM*A (K > 

35  CCMIMJE 
«PkCC=SUM 
Kiltkh 
cMM  ASUR 
IMM  GO.CE’.AC 
Cu  NN- ~  N 
1«  =  AU> 

GC  01  K - i 
C  (K  I  -AAH)  (K  t 
0  1  CCM  INUt 

0  2  Ht,  1  Lots 
00  UC  AS  K= 1  ,N 

COO-AOO  «UOO 
08  CCM  IMJE 
ht 1 CRA 
cl.tRY  5HAP 
OC  50  K= 1 ,N 
AA-AIK) 

AOO=b(K) 

U  (K)  =  AA 
5  0  CCM  1  MJt 
R  i  I  C  Ob 
i.MRY  001 
111  0. 

lie  CO  K-  1  ,N 
Of  CP»AOO*UOU 
m;  clMIMJe 
•m-aCC-DP 
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nt  I  URN 
tMH  XM 1 1 
16 <N>  b2,fcb,6 7 
b2  NN  =  -N 
A  A=  A  111 
Dc  fcb  K=1,NN 
fc  IK>  =  AA 
fe5  CCMlNUE 
66  he  1  CAN 
fc?  00  70  K=1.N 
E  ( K  )  -  A  ( <  > 

7 it  CCN11KUE 
fit TUAN 
fc  N  C 


j 


n  n  n  n 
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FUNCTION  GAMA  (7) 

COMPUTES  GAMMA  FUNCTION  OF  A  REAL  ARGUMENT 
SEE  FEF.4,  ECS  6.1.15,  6.1.16,  AND  6.1.35. 

DATA  Al/- . 5T 46646/ , A  «/• 951 2363/ , A3 /- .6556566/, A4/.4 24 5545/ , 
1  A5/” « 1 0 10  67  8/ 

X-7-1. 

IF1X.LT.0.1  STOF 
GAMMA-1. 

ia  if«x.lt. i.i  r,o  to  ?o 

GAMMA -GAHMA*X 
X-X-l. 

GO  TO  10 

20  GAMMA-GAMMA • «U ( ( A5*X«A41 "X ♦ A3 1  *  X* (2 > «X* Al » *  X* 1 . I 
RFTUPN 
£  NO 


SU»ROUTINE  FRNTIX,N,NPL,A,KHEREI 


C  THIS  SUBROUTINE  PRINTS  OUT  ARRAY-X  OP  LENGTH  N.  THE  FORMAT  IS . 

r 

C  NPL.CT.O  F-FORMAT,  NFL  CECIMAL  FLACFS 

C  NPL.LT.C  E-FORMAT, -NFL  CECINAL  PLACES 

C 

C.  IN  AOniTICN . 

r 

C  A  =  HOLLERITH  NAME  OF  ARRAY  BEING  PRINTED  (6H,  RIGHT  JUST  I  FIE Cl 

f  HHERE  -  HOLLERITH  LABEL  (20H1 

r 

DIMENSION  X ( II »NHERE 12Y 

PRINT  100,A,N,WHEF£,A, (A,K,K=1,91 

T  1  =  1 

12  =  M INC (N, 10  * 

10  NF  = I ABS (NFL  1 

IFINFL.GT.OI  PRINT  101, II, NP, (XIII ,1  =  11,12  I 
IMNFL.LT. 31  PRINT  102, II, NP, 1X111,1  =  11,121 
It=Il*lC 
I2  =  NJNC'1I1*<3,N'> 

IFlH.Lt.NI  GO  TO  10 
RPTURN 

100  FORMAT  I////13H  PRINTOUT  OF  , A6,8H1KI ,K=1,IL. CX2A10// 

1  5H  K,6XAG,.THlNI,qilXA6,3HlN4,Il,lHI  I/I 

101  FORMAT (IS, IX10F12.=| 

102  FORMAT (15, ?X 10E 12. =1 
ENF 


sufrgutine  spctrm ixr , xt  .ntn.nout .alfha.wi 

r 

C  IN  THIS  StfiFOHTINE  WE  COMPUTE  THE  FGKER  SFECTRUM  OF  TEE  COMFLEX  TIKE 

c  scquence  in  the  array-fair  cxr.xit  of  length  kin.  the  power  sffctfck 

C  IS  RLTURNEC  IN  ARRAY  Xfi,  ANC  IT  IS  NOW  OF  LENGTH  NOUT. 

r 

C  THC  SAMPLE  SPACING  OF  THE  FOHER  SPFCTRLM  IS  1/NOLT  OF  THE  REPETITION 
C  FPEGUENCY. 

C 

C  A  COSINE-ON-A-FEOESTAL  WEIGHTING  IS  APFLIED  TO  the  INFUT  SAKFIFS. 

C  ALPHA  IS  THE  PATIO  OF  THE  HEIGHT ING  FUNCTION  AT  THE  ECGE  TC  THE 
C  CENTER.  ALFHA=.0E  FOR  HAMMING  AND  ALPHA=1.8  FCR  UMFCRM  WEIGHTING, 
r  THE  USER  CAN  SUFFLY  HIS  CNN  WEIGHTING  FUNCTION  IN  ARRAT  W  CF  LENGTH 
C  NIN.  RUT  HF  MUST  SFT  ALPHA  TO  A  NEGATIVE  VALUE. 

C 

C  AFRAY  W  IS  A  WORKING  ARRAY  ANC  IT  MUST  RE  DIMENSIONEC  AS  LARCE  AS.... 
C 

C  NIN  IF  NOL'T.EG.E**INTEGER 

C  NTN42*N0UI  IF  NOUT. NE. 2* "INTEGER 

C 

C  THE  Wfir.HTS  ARE  NORMALIZED  SO  THAT  THE  SUM  IS  LNTTY. 

C 

CIKENSTCN  XR(lt,XI(l).N(ll 
DATA  NTN9/0/ .NGUTC/O/ 

IFtNIN.EO. NINO. AND. NGUT .EC.NOUTO  >  GO  TC  1C 
IFIALFFA.GS.0.1  CALL  HEIGHT tW, NIN, ALFHAI 
NINO-NIN 
NGLT  0  =  NGUT 
19  N-NOUT-MN 

IFTN1  <,0,15,12 

12  CALL  XMIT|-N,a.,XRCNTN»l)l 
CALL  XKITC-N,a.  ,XI  (NIN  +  DT 
IF  CALL  FROniMN.H.XR.XPI 
CALL  PROCMMN.H.XT.XIT 
CALL  FFT2 (>R,XI,NCUT,-ll 
CALL  PCWRINOLT.XP.XI.XRT 
RCTUPN 

<,0  FFINT  ioo.nin.ncut 

sror 

100  roPMATI/lX<,9H»4»EFH0R  IN  SPECTF.K,  NOUT  IS  SMALLER  THAN  MN  •*•/ 

1  /lnxRHNIN'IR , 10XSHNOUT  =  T  M 


[-  -AD-A097  765  MARK  RESOURCES  INC  MARINA  DEL  REY  CA  F/G  17/9 

SIMULATING  RADAR  SIGNALS  FOR  DETECTION  PERFORMANCE  EVALUATION. (U) 

FEB  01  R  L  MITCHELL  F19628-79-C-0012 

RADC-TR-BO-410  nl 


UNCLASSIFIED 
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SUFROUTINE  HEIGHT (W,N, ALPHA! 

THIS  SUBROUTINE  COMPUTES  CGSIKE-ON-A-PEOESTAL  HFIGHTS  IN  AKRAT  H  CF 
LENGTH  N  «N  IS  ALSO  THE  NUNBFR  OF  PULSES*.  ALPHA  IS  THE  RAT1C  OF  THE 
WEIGHTING  FUNC T IGA  AT  THE  EOGE  TG  THE  CENTER.  ALPHA-. 0«  FCR  HANNING 
ANO  1.0  FOR  UNIFORH  WEIGHTING.  THE  WEIGHTS  ARE  NORN ALI7ED  SC  THE  SUN 
IS  UNITY. 

IF  THE  WEIGHTS  ARE  TC  FE  LSEO  AS  PART  CF  A  LONGFR  ARRAY  OF  LENGTH  NR 
(WITH  7ERO  FILL!.  THE  CALLING  SECULNCE  IS . 

CALL  XMITI-NR  ,0. ,WI 
CALL  WEIGHTSTW.N, ALPHA* 

TO  CENTER  THE  WEIGHTS  AT  THE  FIRST  SAHPLE  CF  THE  W-ARRAV  fSESI  DONE 
WHEN  N  IS  C001,  FOLLOW  THE  ABOVE  SEQUENCE  WITH . 

CALL  SNAFT»W,NP,-N/2» 

OINENSION  kill 

OATA  TWOFI/6. 2631553/ 

»=»l.*ALFHAI/2. 
n  =  ( 1 . -  ALPHA  I /2 • 

r.NrtN»l>/2. 

XN  =  N 

CO  20  K  =  1,N 

W<M  =  A*fl*COS<TWGFI»<K-CN»/XN! 

20  CONTINUE 

HN0Rh=l./7SUK»N,H> 

CALL  FCCOt-N.WKCRH.W.WJ 
HFTUFN 

f.  nc 


oononnoo 
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FUNCTION  CeXJA.JU 

COMPUTES  OECI  PEL  TALVfcS.  THERE  APE  TMC  ENTRIES . 

0=DFIA»  d=ib.*alogio iay 

CALL  OON  IA.M  •  A  IKY  =CB  (A  IK  1 1 «  K=1,N 
THE  OUTPUT  IS  TRUNCATEO  TO  THE  INTERVAL  -RS,  C». 
CUNFNSTON  Mil 

nE7|7»=10.*ALOCia t AMAX1 « AK1N1 17. «AJES .2»  *1 «258SE-10n 
fNTRV  OB 
r.BX:nB7IA(11l 
FfTURN 
ENTRY  CBN 
CO  70  1=1. N 
A 1 1 1 =DE7 ( A  fill 
70  CONTINUE 
RETURN 
FNC 


y  i 


SUBROUTINE  DSTGEN tA.P.NNI 

COMPUTER  CUMULATIVE  OISTRIEUTICN  FUNCTION  CF  SAMPLE  CAT* 

TO  INITIALIZE . 

CALL  OSTIAL  1X1, XINC.Kt  XI  -  FIRST  VALLE 

XINC  =  VALUE  IKCFEHFNT 

N  =  NUKPfR  CF  INCF'KENTS 


IOR  EACH  CATA  POTNT . 

CALL  OSTPNTIX.NI 


X  -  DATA  PCINT  VAl  JE 
H  =  HEIGHT  USUALLY  *1.1 


TO  COMPUTE  THE  GISTPIflLTION  FUI  'ION  (AFTER  LAST  CALL  TC  CSTFNTI.. 

CALL  CSTFUMPT  F  =  OISTRIEUTICN  FUNCTION 

THE  CIETRir.UTICN  FUNCTION  IN  ARRAY  F  15 . 

Mil  =  FRCntOATA.LT. XII 
PI2)  =  FPCniTiATA.LT.Xl  +  XlNCI 


C 

C  FfNI  =  FRCetOATA.LT.Xl*tN-ll»XINCI 

C 

C  ARRAY  F  IS  riMENSICNEO  FOR  A  MAX  VALLE  OF  N=201 
C 

OTMENSTON  A  til  ,F (E02I 
ENTRY  rSTINL 
N  =  NN 

IPtN.GT.EDU  STOP  22 

NS  =  0 

NPi -N*l 

X1-AUI 

v  T  NC  -  P. 

CALL  XPITt-NPl,0.,FI 
RETURN 

ENTRY  r.STFUN 
TO  10  1=2, NPI 
rtLi-rti-ii  +  FUY 
10  CONTINUE 
HO  20  L  =  1  ♦  N 
F  tL  1  -  H (NF 1  » -F tLT 
A(LI  =  1.-P ILI /FLOAT  <NSI 
20  CONTINUE 

C  PRINT  100  «P  tNRll 

iau  roFMATtrift.fei 
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RETOOK 

entry  ostpnt 

NS=NSM 

L=IM1»-X1>/XIKC*Z. 

l  =  M«X0  (L « 1 1 

L-=HlM>  TL.NF1T 

OU.»=ETL»*F 

RETURN 

FNC 


nomnnom 
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SUBROUTINE  FFT2tAR,AI,N.ISCN) 

THIS  SUBROUTINE  COMPUTES  THE  FFT  IN  THE  ARRAY-PAIR  IAR.AII  CF  LENGTH  K 
ANO  RETURNS  THE  RESULT  IN  THE  SAME  ARRAY  PAIR. 

I SGN  =  SIGN  Cr  FHASE  ARGUMENT  IK  FPT 

ARRAY  M  IN  COMMON  MUST  BE  OIHENSIONEC  AT  LEAST  AS  LARGE  AS  ?»N  IF 
N.Nt.2»* INTEGER. 

OIMENSION  AR  (II . AT  (1 ) 

COMMON  /TEMP/  Hill 
1  =  1 

IFIISGN.GT.OJ  1=3 

CALL  FCUfiTIAR.AI.Ktl.IfltHdl  .HIN41M 

RETURN 

ENC 


(JUUOtUU  C-O  U  U  C.  w  C. 
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SUBROUTINE  FCURTIOATAR.OATAI,NN,KC!M,IFRWr,irFlX,WCRK*,NCB KI» 
DIMENSION  E A  TAR «1 I,  CATA I (II .NNI11 ,MOh*R < 1) , MOPKI tl» ,1  FACT  120 > 
f 

C  THE  COOltY-TUKEY  FAST  FOURIER  TRANSFORM  IK  USASI  BASIC  FCRTRAN 

C 

f  EVALUATES  COMPLEX  FOURIER  SERIES  F CR  CCKPLEX  CR  REAL  FUNCTIONS. 

C  THAT  IS,  IT  COMPUTES 

C  FTRAMJl.J2,...»»SUMICATAm,T2,...»«tll«*tll-lT*IJl-n 

f  *N2**tI2-ll*«J2-l»*...»» 

C  WHERE  Nl=EXPf-2*FI*SCRTt-lT/KNflll,  N2=EXP  t-2*  FI  *SCRTI-i  t/NN'21 1 , 

C  ETC.  AND  II  ANO  J1  RUN  FROM  1  TO  KK*1».  I?  ANO  J2  RUN  FRCK  1  TO 

f  NN(2T,  FTC.  THERE  IS  KO  LIMIT  OK  THE  DIMENSIONALITY  (NUMBER  OF 

C  SUBSCRIPTS*  OF  THE  ARRAY  OF  DATA.  THE  PROGRAM  MILL  PERFCRM 

C  A  THREE-DIMENSIONAL  FOURIER  TRANSFORM  AS  EASILY  AS  A  ONE-GIMEN- 

C  SIONAL  ONE.  TKO  IK  A  RfiOFCRT lONATEtY  GREATER  TIME.  AN  INVERSE 

C  TRANSFORM  CAN  BE  PERFORMED,  IN  MUCH  THE  SICE  IN  TEE  EXPONENTIALS 

C  IS  ♦,  INSTEAD  CF  IF  AK  INVERSE  TRANSFORM  IS  FESFORMEC  LPCN 

r  AN  ARRAY  OF  TRANSFCRMEC  PATA,  THE  ORIGINAL  DATA  MILL  PEAFFE**. 

C  MULTIPLIED  BY  KM  111  *NN  (  Zl  »  , .  THE  ARRAY  OF  INPUT  TATA  MAY  BE 

r  PEAL  OF  COMPLEX,  AT  THE  FECGRAHMERS  OPTION,  WITH  A  SAVING  CP 

C  ABOUT  THIRTY  PER  CENT  IN  RUNNING  TIME  FCR  REAL  CVEF  COMPLE  >• 

C  (FOP  FASTEST  TRANSFORM  OF  REAL  DATA,  NN tit  SHOULC  BE  EVEN. I 

C.  THE  TRANSFCRM  VALUES  ARE  ALWAYS  COMPLEX ,  AND  ARE  PETURNEC  IN  THE 

C  ORIGINAL  ARRAY  or  CATA,  REPLACING  THE  INPUT  DATA.  THE  LENGTH 

C  OF  EACH  DIMENSION  GF  THE  CAT t  ARRAY  MAY  PE  ANY  INTEGER.  THE 

C  PROGRAH  RUNS  FASTER  CN  COM,FOSl*£  INTEGERS  THAN  CN  PRIMES.  AND  IS 

C  PARTICULARLY  FAST  ON  NUMBERS  RICH  IN  FACTORS  CF  TKC. 

C 

C  TIMING  IS  IN  FACT  GIVEN  BY  THF  FCLL0WTN6  FORMULA.  LET  NTCT  «E  TME 

C  TOTAL  NUMBER  OF  RCINTS  IREAL  OP  COMPLEXt  IN  THE  DATA  ARRAY,  THAT 

C  IS,  NTCT=NN«l»*NNt2l»...  CECOMPCSS  NTCT  INTC  ITS  FRIME  FACTORS, 

C  SUCH  AS  2**K2  •  3*«K3  •  S**K5  *  ...  LET  SUM?  RE  THE  SUM  CF  All 

C  T HR  FACTORS  OF  TWC  IN  NTOT «  THAT  IS,  SUM?  *  2**2.  LET  SIM.F  ef 

C  THF  SUM  OF  ALL  OTHER  FACTORS  OF  NTCT,  THAT  IS,  SUM*  s  3*K3*5*K5«.. 

r  THE  TIME  TAKEN  BY  A  MUL T I G I MENSI CNAL  TRANSFORM  GK  THESE  NTCT  CATA 

IS  T  :  TM  Tl*NTOT  «  T2*NTOT*SUM2  ♦  TT*NTCT *SUMF .  FOR  THE  PAP- 
TTCULAR  IMPLEMENTATION  FORTPAN  32  CN  THE  CCC  7300  f FLOATING  POINT 
AGO  TTMC  =  SIX  M.IfROSECONCSI  , 

T  -  3000  ♦  600*NTC,T  ♦  50*N10T»SUM2  ♦  1 70  *NTO  T*  JUKF  MICROSFCCNOS 
ON  COIFLEX  nATA. 

IMPLEMENTATION  CF  the  CEFIMTIOK  RY  SUMMATION  Will  PUN  IN  A  T IMF 
POCF&RTIONAL  TC  NTCT**2.  FOR  HIGHLY  CCNPCSITE  NTCT,  THE  SAVINGS 
OFFEREE  BY  CCCLFY-TUKFY  CAN  BF  OFANATIC.  A  MATRIX  100  EY  100  WTll 
r-F  TFANSFPRHO  IN  TIME  PRCFORT  TONAL  TC  100  00  *  C  2*  Z*  c  *Z*  St*  ♦  S*S>  = 
2B0.00C  (ASSUMING  T2  ANO  T3  TO  BE  ROUGHLY  CCMPARAFVct  VERSUS 
lCC'Oi  **?  -  lOO.dOO.OOC  FOR  THF  STR  A I GHTFOR  WAFB  TECHNIQUE. 

T Hr  CCCLEY-TUKET  ALGCRITHN  PLACES  TWC  RESTRICTICNS  UPON  THE 
NATURE  OF  THF  OATA  BEYOND  THE  USUAL  RESTRICTION  THAT 


0100 
:ioi 
0102 
0103 
010* 
OIOS 
010G 
0107 
010* 
0109 
0110 
Clll 
0112 
0113 
01  IN 
011S 
C116 
Oil? 

one 

0119 

5120 

5121 
0122 
0123 
012* 
C12S 
0126 
012T 
012* 
0129 
0130 
0131 
0132 
0133 
013* 
011S 
0136 
013T 
"138 
5139 
01*0 
01*1 
01*2 
01*3 
J!  ** 

n*s 
?1*6 
01  *7 
01*8 
01  *9 


% 


C  the  data  FRO*  ONE  CYCLE  of  A  PfRIOCIC  PUNCTICN.  they  ARE--  C150 

F  1.  -HE  NUMBER  OF  INPUT  OATA  AND  TEE  NUMBER  C*  TRANSFORM  VALUES  0151 

C  MUST  »E  THE  SAME.  C152 

C  2.  CONSIDERING  T*IE  OATA  TO  BE  IN  THE  TIME  OCMAIN,  0153 

f  THEY  MUST  be  ECUI-SPACEC  AT  INTERVALS  CF  OT.  FURTHER.  THE  TRANS-  0154 

C  FORM  VALUES.  CCNSICEREC  TC  RE  IN  FREOtEKCY  SPACE.  RILL  BE  EOUI-  0155 

C  SPACIO  FROM  0  TO  2*PI* (NN !I1-11/(NN!I1*0TT  AT  INTERVALS  CF  015ft 

C  2*PI/  INM  IT  •OTV  FCR  EACH  CIHENSICN  OF  LENGTH  KMII,  OF  CCtfiSE,  C1S7 

C  OT  NtEr  NOT  BE  THE  SAME  FOR  EVERY  DIMENSION.  C158 

r  0159 

C  THE  CALLING  SECUENCE  IS--  0160 

C  CALL  FCURTCOATAR, OATAI, NN.NOIM.IFRWD.ICPLX.WORKR.WCRKIT  0161 

r,  C162 

f  DATA*  ANC  DATAI  API  THE  ARRAYS  USED  TC  HOI D  THE  REAL  ANC  ’MAGINARY  0163 

C  PARTS  OF  THE  INPUT  OATA  ON  INPUT  AND  THE  TRANSFORM  VALUES  ON  0164 

C  OUTPUT.  THEY  ARE  FLOATING  POINT  ARRAYS,  MUL  II C I  MENS I  GOAL  KITH  0165 

C  IDENTICAL  C I MENSI CNALI I V  ANO  EXTENT.  THE  EXTENT  CF  EACH  OIKENSICN  0166 

C  IS  GIVEN  IN  THE  INTEGER  ARRAY  NN,  CF  LENGTH  NCIN.  THAT  IS.  0167 

f  N01M  IS  THE  01 PENSION ALITY  OF  THE  ARRAYS  DATAR  AND  DATAI.  0168 

C  IFPWO  IS  AN  INTEGER  USED  TC  INGI CATE  THE  OIFECTICN  OF  THE  FOURIER  0169 

C  TPANSF CRH.  IT  IS  NCN-7ERC  TC  INCICATE  A  FCRWARD  TRANSFORM  OlTC 

C  l£X  °ONENT I AL  SIGN  IS  -I  ANO  7ERO  TC  INCICATE  AN  INVERSE  TRANSFORM  0171  1 

f  (SIGN  IS  ♦>.  ICFLX  IS  AN  IN TfcGFR  TO  INCICATE  WHETHER  THE  CATA  0172 

C  ARE  PEAL  OR  COMPLEX.  IT  IS  NON-7EPO  FCR  CCMFLEX ,  7ERO  FCR  REAL.  0173 

r.  IF  IT  IS  7FRO  TREALI  THE  CONTENTS  CF  ARRAY  OATAI  HILL  EE  ASSUKEC  0174 

C  TO  Bl  7EKO,  ANO  NEED  NCT  BE  EXPLICITLY  SET  TC  7ERC.  AS  EXPLAINEC  C1T5 

C  AEOVE,  THE  TRANSFCRM  PESULTS  ARE  ALWAYS  COMPLEX  ANC  ARE  STCREC  0176 

f  IN  OATAF  ANC  OATAI  CN  RETUPN.  WCRK®  ANC  WCPKI  ARE  ARRAYS  USEC  0177 

f  FON  WORKING  ST GR AGE .  THEY  ARE  NCT  NECESSARY  IF  ALL  THE  CTMENSICNS  0178 

f  OF  THE  OATA  ARE  POWERS  OF  TWO.  IN  THIS  CASE,  THE  ARRAYS  HAY  BE  0179 

C  REPLArtf.  BY  Tnr  NUMBER  0  IN  THE  CALLING  SECUENCE.  THCS.  LSE  C F  0180 

r  POWERS  cr  TWO  CAN  free  a  GCOC  DEAL  CF  STCRAGE.  IF  ANY  DIMENSION  0181 

C  IS  NOT  A  POWER  OF  TWO,  THESE  ARRAYS  MUST  BE  SUPFLIEO.  TH;  Y  APE  0182 

r  FLOATING  PCINT,  one  G.IMENSIONAL  CF  LENGTH  EOUAL  TC  THE  LfJCFST  0163 

C  ACFAV  r I  HE  NS  TON,  THAT  IS,  10  THE  LARGEST  VALUE  OF  NNf II .  0164 

C  WCRKfi  ANO  NCRKI,  IF  SUFFLTEO,  MUST  NCT  BE  THE  SAME  ARRAYS  AS  CATAR  "185 

r  OR  OATAI.  ALL  SUBSCRIPTS  OP  ALL  ARRAYS  BEGIN  AT  1.  0186 

C  0187 

C  EXAMPLE  1.  THFEr-SIMf NSICHAL  FORWARC  FOUR IFP  TRANSFORM  C  A  3188 

C  r.OHPlr>  ARRAY  G I  HENS  I ONEO  100  EY  If  RY  13.  0185 

f  OIMFNSTON  rATARUOO. 16.131, nATAKlOO. 16. 13T.WORKRI100).NCRKI(100I  3190 

G  OIMrNSION  NN(3»  0151 

r  NNIt I  SCO  0152 

G  NNI2I-16  0153 

G  NN! ’ 1 -13  0154 

C  r.AU  rGURT  (TATAR, CATAJ  ,NN,3,  1 , 1  .  NOR  KP ,  KC  RK  IT  0155 

r  0156 

G  CXAM1L  ?.  <'NL  -  0  I  “F  NS  1  ON  AL  FORM  ART  TRANSFORM  CF  A  REAL  ARRAY  OF  0157 

G  IEN<  TH  64.  0158 

r  r  I M*  N’T  ON  r  A  rAM64»  ,r  ATAT  !64l  0159 


% 
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f  CALL  rcuKT  « 0  AT  AR ,  OAT  A I  ,  64,1.  1,D,  C.9Y  0?00 

r  C20i 

C  there  ARE  NO  FRROP  MESSAGES  OR  ERKLR  HALTS  IN  THIS  PROGRAM.  THE  C202 

C  FBOGRAH  RETURNS  I MMEOIATf LY  IP  NOIN  CR  ANY  NNIII  IS  LESS  THAN  ONE.  5203 

r  9204 

C  The  SINE  ANO  COSINE  VALUES  REQUIRE!  FCK  THE  TRANSFORM  ARE  0205 

C  GENERATE 0  RECURSIVELY.  IF  OOURLE  PRECISION  IS  A  VAILAELE «  IT  IS  C206 

C  STRONGLY  UR6E0  THAT  THE  FOLLOWING  VARIABLES  «E  SC  DECLARED  TO  5207 

C  RECUCF  ACCUMULATION  CF  ROUNDOFF  ERROR--  0208 

C  DOUBLE  PRECISION  TNOP I , THE 1A .WSTFR , NSTFI ,WNI HP ,KHINI, MR, VI ,NTf HP  0209 

r  •  ,THETH,HHSTR,HMSTI,TMOHR,SR,SI,CLOSfi,CLOSI,STHFR,STKFI  0210 

C  IN  ADDITION,  THOPI  SHOULD  EE  ASSIGNFC  A  SUFFICIENTLY  PRECISE  0211 

C  VALUE  ANO  THE  VARIOUS  CALLS  TC  THE  FUNCTIONS  CCS  ANO  SIN  0212 

C  SHOULD  BE  CHANGED  TO  DCOS  AND  DSIN.  0213 

C  C214 

C  PROGRAM  BY  NCRMAN  F.RENNER  FROM  THE  BASIC  ALGOR  ITH'H  BY  CHARLES  C215  ! 

C  RADER  TOOT H  OF  NIT  LINCOLN  L ABORATCR Yl •  HAY  1967.  THE  IDEA  5216 

C  FOR  THE  BIT  REVERSAL  MAS  SUGGESTED  RY  RALPH  ALTER  TALSC  HIT  LL».  0217 

C  ADAPTEC  FRCH  the  MORN  OE  JANES  M.  COOLEY  ANO  JOHN  Hi.  TUKEY .  C210 

C  AN  ALGGFITHH  PC®  THE  MACHINE  CALCULATION  OF  COMPLEX  FOURIER  0219 

f  SERIFS,  MATH.  COMPUT.  19,  90  AAPRIL  1=651,  2=7-301.  C220 

C  0221 

IFINTTM-H =20,1,1  C222 

1  HTOT  =  1  0223 

DC  2  T OI M= 1 , NO IM  0224 

2  NTOT=:NTGT*KNLiriM>  0225 

THOPI- G.2831A53C 7  0226 

9  22T 

MAIN  LCCF  FO=  F ACH  CIMENSICN  0228 

0229 

NF1  =  1  32  30 

CO  913  miH^l.NOIH  0231 

N=NNIiriK»  3232 

NF2rNFi*N  !<233 

IFCN-11920,900,5  0234 

0235 

IS  N  A  POWER  OF  TMC  ANC  IF  NOT,  NH A T  ARE  ITS  FACTORS  0236 

p  2  37 

H-N  0238 

NTMONF1  0239 

TF=1  0240 

IDTV  2  0241 

10  lOUPT-M/iniV  n?l,2 

TFIM-M-IOTVOTCUOT  0243 

lFIICtlUT-Tnvi50.ll.il  7244 

11  IFITPFMI20. 12,20  n?i,5 

12  NTVC  NTHCmMNO  n2*f, 

TFACI * IT1 : IOIV  0247 

IF-TFM  0248 

M  =  IGUC  T  C249 


20 


9H 


r,o  tc  n  o?50 

iciv^s 

IKCN?=IF  3252 

30  I  GUOT=H/I 01 V  0253 

IREM=H-IOIV*IOUOT  0254 

IFl IOUCT-ICI VI 60. 31. 31  025$ 

31  IFCIREMAO. 12.40  3256 

32  IF ACT IIFI=IDIV  5257 

IF=IF»1  025* 

H=TOUOT  5259 

GO  TG  30  0260 

90  IOIV=irlV»2  5261 

GO  TO  ?C  0?6? 

50  INCN2- IF  0263 

IF1IREMI60.51.60  0264 

51  NTWO=NTKG*NT«G  0265 

GO  TO  70  5266  , 

60  I F  ACT  t IFl  =  M  0267  | 

70  NON2F - NP2/NTW0  026# 

C  0269  I 

r  SEPARATE  FCUS  CASES —  5270  ‘ 

r  1.  fOMPLEX  TRANSFCPM  0271 

C  2.  REAL  TRANSFORM  FOR  THE  2ND,  ?RC,  ETC.  riMENSION.  EETHGC —  0272 

c  TRANSFORM  HALF  THE  CATA,  SUPFLYING  THE  GTHE F  HALF  EY  CCK-  0273 

C  JUGATE  SYMMETRY.  0271. 

C  7.  crflL  TRANSFORM  FCR  THE  IS."  DIMENSION,  N  CCC.  HETHCT--  C275 

r  sn  the  IMAGINARY  PARTS  TO  7ERC.  0276 

C  A.  r F AL  TRANSFORM  FGR  THE  1ST  OIHENSION,  N  EVFN.  METFCt--  3277 

r  TRANSFORM  A  COMPLEX  ARRAY  CF  LENGTH  N/2  WHOSE  REAL  FARTS  027# 

r  ARE  THE  EVFN  NUH8EREC  RFAL  VALUES  AND  WHOSE  IMAGINARY  FAFTS  0279 

C  ARE  THE  OOO-NUMFEREC  REAL  VALUES.  UNSCPAHPLf  ANC  SLFPLT  0280 

r  THE  SECOND  HALF  EY  CONJUGATE.  SYMMETPY.  0281 

C  0282 

TCASF-1  0283 

I FMT  N : 1  0284 

IF1TrTixilOO.71.100  3285 

71  ICASr-2  02*6 

TFUriE-ll 7E,7?. ICO  0287 

73  IF  A  S  E  r  3  0288 

TF1NTWC-NF1 HOC. 1C0.7T  0289 

73  I CA  tt - U  52S0 

IrHTN-2  0291 

HTWO-MWO/2  0292 

N-N/2  0257 

NP-  M.-/7  9294 

NTOTNTGT/2  0295 

T  =  1  0296 

no  ca  j-i  .mot  0297 

rATAP  rjl-OATARtn  0298 

OAT  AT  ulrnATAR  TIM  )  0299 


I  -  I *2 


')') 


no 

r 

r 

c 

r 

100 

101 


121 


125 
13C 
14  0 
11.1 


153 

r 

c 

r 

200 


210 

220 

24C 

250 

25'. 

21*0 


SHUFFLF  DATA  Ot  BIT  fif VEkSAL ,  SIFCE  Nrj***.  AS  THE  SHUFFLING 
CAN  BE  CONE  FT  SIMPLE  INTERCHANGE,  NO  KCRKING  ARRAY  IS  NEEOEB 

IF 1NON2F-11 101,101,200 

NF2HF  =  NF2/2 

J=1 

OG  lr*0  T2-1.NP2.NF1 

IF  1.1-12)121, 130, 130 

I1HA*-I2*NF1-1 

GG  125  Hr  12  ,I1KAX 

00  12r*  I3rH  ,NT0T,NO2 

J'=.J*I3-I2 

TENf'FrrATAMIS) 

T  CHI'  IrfATATII.T) 

[3ATAMI3)rCATARTJ?T 

DATAIII3)=CATAIIJ3* 

0ATARIJ3)rTEHPF 
OATAHJ3I -TERPT 
M=K02NF 

IF(J-“)  150, 150,11.1 
J  =  J-  F 

HrH/2 

IFTM-Nr- 11160, 140, 140 
J  =  J*  F 
GO  TC  300 

SHUFFLE  DATA  BY  DIGIT  REVERSAL  FCR  GENERAL  N 

00  273  Ilrl.NPl 
00  270  I?rTl  ,NTOT  ,NPZ 
J=  I  3 

GG  2F0  T~ 1 , N 
IF(ICArE-3)210,220,210 
RORKR IIlrOATARTJI 
WORKI(I)rCATAUJ) 

GO  TG  240 
UORKMI)  "DAT  ART  J) 

HGPKI  (T1=0. 

I F  F2  =  NF2 
IFrIFPTN 

I  FF 1 r I rF2  7IFACTITF) 

Jr Jtifri 

IF! J-l »-IFF2(?GU ,255,?rq 
JrJ-IFF? 

[FI  ?rlf  ri 

C  F  -- 1  F  ♦  1 

!  r  t  T »  I2-NP1 t  750,  ?F.0,2r0 
.OM  IMIE 


9390 

C301 

0302 

0303 

0304 

0305 

0306 

0307 

0300 

C309 

0310 

C311 

3312 

0313 

0314 

0315 

0316  ! 

0317  ; 

C318  ! 

0315  1 

0320  ! 

0321 

0  322 

0323 

5324 

0325 

0326 

032-F 

0328 

0  320 
0  3  30 
0331 
0332 
0333 
0334 
5335 
C  336 
0337 
0  138 
C  339 
0  3  40 
C  341 
0  342 

9343 

0  344 
0345 
0  346 
0  347 
0  3  45 
0  349 


100 


I27AX=I3*NF2-NF1 

0350 

1=1 

0351 

00  270  12  =  13  ,I?MAX,NP1 

C  3  5? 

0*TAfi»121=HCKICMIl 

0353 

0ATAI(I21=WC9KI1II 

0  3  51, 

270 

1  =  1*1 

0  3  55 

r 

0356 

c 

SPECIAL  CASE—  W=1 

C357 

c 

0358 

3QC 

1 1PNG=6F1 

0359 

GO  TO(302,»01,302,3C2>  .ICASE 

0  =  60 

301 

llF6G=6F0*ll+NFREV/2) 

C  361 

302 

IFINTWC-NP1 1600,603, 303 

0362 

303 

CO  430  11=1,  I1FNG 

0363 

I  )*IN=NP1  ♦  1 1 

0364 

1  STEF  =  2  *NP1 

0365 

GO  TO  330 

0366 

310 

J  =  T1 

0  36? 

CO  320  1=  16  IF, MOT  ,  ISTE F 

0368 

T  EMPF  =  C AT  AR 1 1 ) 

0369 

T  CMP  I  =  r  AT  A  I  in 

0  370 

r ATARI 1)=0AT API J» -TEHPC 

0371 

r  ATfll  (  H-PATAI  (J1 -TE7FI 

C  372 

1  AT  AT  1  J1=0AT  AR  131  *Tt  hPR 

037  = 

DATA  11 J»=OA  T»I 1 J»  *  TEMPI 

0374 

320 

J=J*ISTtP 

0375 

IHTN=TPIN* IMIK-ll 

0376 

ISTEF  ISTEFMSIEF 

3377 

330 

IF (1ST f  F -NT KOI  310,310 .331 

037? 

r 

0  3  7° 

r. 

SPECIAL  CASE—  H  =  -SORri-ll 

C  3  80 

r. 

0381 

331 

I HIN -  3  *Npl <11 

03*2 

TSIEf -4*NP1 

0  38  = 

r.O  TO  4?0 

0384 

40C 

IMMIN-ISTEF/2 

C385 

TO  410  I=I('IF,KTOT,ISTEF 

0  386 

IF  (I)  FWC143  1  ,41.?  ,  401 

0387 

<•01 

t  f MPF-pAT  A  T 1 II 

0388 

T  EF-r  I  -CATARII1 

0389 

r,c  tc  40  3 

0390 

10  2 

T.KPF  -CATAIin 

0  391 

TEMPI  TATAR! I) 

:  =  52 

40  3 

0  AT  A  k  1  I  1  -r.  ATAO  «J>  -U  mpf 

0393 

'iATAI  1  T  >  =DA  TAI 1J) -TEM  I 

a  394 

CAT  AF  1 J) -DATAF1 J)  ♦TEH'S 

'  39S 

0  AT A I  1 Jl-TATAIHJ) *T£  EPl 

0  3  96 

410 

J=J*IST£P 

0397 

I  MIN  ■  IMNMMIN-Tl 

0398 

ISTf.f  IMF  T  ,  ISTtF 

0399 

101 


52C 

MG 

r 

r 


r 


501 

50? 

CBC 


rlC 


i 

[  *?0 


5?5 

mo 


550 

r 

C 

r 

r 

( no 

coi 


IMISTEF-NTWOI50  0.500.530 
CONTIM’E 

FAIN  tOOF  FOF  FACTGRS  Of  TWO.  W  =  E  X  P  1-2»PI  *SOP  T  1  -1 »  •H/t'KtX  I 

THETA--THCPI/8. 

WSTFC-o. 

wsTri=-i. 

IMTFFVGiSO?  .501,50? 

THETA=-THFTA 
WSTf  1=1. 

HKAX=8«Nf  1 
GO  TO  550 
WMlNS=COS(THFTAT 
WHINI  =  SINUHETA» 

WR  =  WM INF 

WIrWWIM 

MHlN="NAX/2.NFl 

ESTFF  =  NF1  *KP1 

no  53?  N=M*IN,KMAX .NSTEP 

GO  5?*  T  1  =  1 ,  II PNG 

ISTEFOKAX 

IElN -K«T1 

J=IHIN-ISTEF/2 

rn  SEC  I'lMN.NTGI.tSUF 

TEMPFMATAE  II)  •WR-F.AT  AI  tl  I  »W  I 

TFEFT^CATAR* I) »WI ♦GATAl (I) *WF 

GAT A=  n>=CATAR ( Jl -TEFFF 

GAT  AHH  =  0ATA1  TJT -TFEFI 

GATAP  <  Jl  =  QATAR  (.11  4TEKFF 

GAT A  I ( J)  =  GATAI  (J1 ♦ TEEP I 

J=J+TSTEF 

INTN  =  IEIN*IMN-I1 

15TEf=lSTEr*IST£F 

If  lISTfF-NTwrj1r10»El0,5?5 

CONTINUE 

WTf Ml  - WK»WSTCI 

Hp-H»*l.<-rPF.WI»wc,1Fj 

h  I  t W  I  • 5  M  F  f  ♦  h  T  f  E  r 
WSTPF  NMINF 

w  si  p  r  wf.iM 

THfTt-TNtTA/?. 

MEAX  : E. AX  *EEAX 
IF  TMEAX-NTWOISOO  ,r>  JO  ,f,00 

MAIN  l OOF  F CP  FACTORS  NOT  EOUAL  TO  Two. 

WU  Xf  I  -?*D  I  •  SOFT  t-) ) • ( J2-I31 /If F2» 

If  INCN?F-lt  X  00,  TOO. GO  1 
IfFlrNTWO 


"500 
0531 
050? 
0  50  3 
0505 
3505 
0506 
0507 

0508 
0509 
0510 
C  5  1 1 
051? 
0513 
C515 
0515 
0516 
0517 
0518 
3519 
0520 
3521 
352? 
3523 
0525 
3525 
05  26 
0527 
0528 
f  529 
3530 
3  5  31 
C53? 
0533 
0535 
”535 
0536 
0  5  37 
0538 
''539 
C  5  50 
0541 
0552 
3  55? 
0  5  55 
0555 
3556 

0557 

0558 

3559 


j  y  -  - 
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IF=IACK2 

61J  lrr?=TFACT CIFIMFF* 

1HCTA--TW0F I/FtCAT II  FACT T1F)» 

I F  !  X  F  F  WCT  612  <611*612 

611  THFTA--TMCTA 

612  THFTF  THFTA/FLCATlIFFl/hFH 
HST<*P"C0S(TH”TA> 

hstfi  =  simtheta» 

KMSTF-rCSC  METFI 
HHIT1 -STNUHETMt 
WMIKF-1. 

WHINI-C. 

OP  66C  Jl=l.IFFl,Krj 
UMAX-  J1M1RKG-1 
CO  66 J  11  =  J1  .UMAX 
CC  6*3  11=1 1 ,NTCT  ,KP2 
T-l 

V.R-HMIKF 

KMWK1M 

j?ha»=  1.1‘Ifc?- jf  ri 

00  FAT  .12-1  T,  J2MAX.1FF1 
Twf  ViF  V.S4HR 
JN1N= 17 

JJF'  A  >  =  J2 ’NTJ-IFF? 

CO  630  J3  =  .I2  ,  J’MAX  .IFF? 

J-JKIN«irF2-TFfl 

SR=CATAFtJ> 

SI-PATtl ( J) 

mrsF  =  c. 

OLfGI-O. 

J=j-TFF1 
F?f>  fITMPp  =  5"F 
STMPT  -FT 

6 F rTWOKP^ne-CLCSR^rATAF <J1 
C I ' TWOKF*SI -CLCST  *CATAI < Jl 

PL  GSM  c  TMPC 
0  L  C  C 1 “ 1 Tpci 
J  =  J-l"  1 

II  I  J>. IP  TM  6?  1,621,620 
<  2i  wrcKn'n~wR*sF-wi*3i-oin5;R*nATAP(j» 
WOCKl  U>'WI»SR*«F’'-I.0L0S1  tOATAI  TJ> 
JMTN'.lf  IN*  If  F? 

MI  1  =  1*1 

RT  (  MF  v  R  ♦  w  3  T  F I 
WK  Wf*FSTrt-wl*WMrT 
•  AC  RT  H  I'VSTFF  *  R 1  (  Mr 
T-l 

<  n  f,‘  ii  J2-  I  I  ,irri 

MM  A  X  .l2*Nf,’-T(  e? 

NO  I,'  C  .1*  .1?  ,  J=MAX  .Iff? 


0450 

0451 

0  4  52 

3  453 

1454 

0455 

0456 

0457 

C45« 

3454 

3460 

C461 

C  4  62 

0463 

0464 

0465 

3466  | 

34  6  7  ; 

0468  j 

3469  | 

C 4  70 

34  71 

0  4  72 

3473 

0474 

3475 

0  476 

0477 

34  7? 

C479 

0  4  50 

04  «1 

0482 

0483 

"4  84 
0485 
0486 
0  4  8* 
C4B8 
0485 
A  4  50 
0451 
0452 
3 '•  5  5 
=  1.54 
3  4  5R 
0  4  56 
3.97 
r498 
0<>55 


CATAR(J3I -MORKStH  ?5oo 

PATAI  IJ3l-KCRKItII  0501 

650  1-1*1  C502 

HTEMPsKHINC’tiKSTI  C5Q3 

MMINR  =  WKINR*NMSTP-HHINI*WI>STI  3504 

660  NMIN1-MKINI»WMSTR*WTEMF  0505 

IP  =  IF  +  1  0506 

IFFl^IFr?  C507 

IFJTFF1-NP2I610,700,700  0506 

f  0509 

C  COMPLETE  A  REAL  TFANSFGRM  IN  THE  1ST  DIMENSION,  K  EVEN,  FV  COK-  0510 

C  JUGATE  SYMMETRIES.  0511 

r  0512 

70?  GO  TC  (50 0,600,90«),70ll,IC ASF  C513 

701  NMALF^N  0514 

N-N*N  0515  1 

THETA--TI»OFI/FLOAT(NI  0516 

IF (TFRKG170 3  ,70?, 703  0517 

702  T*ETA=-TFETA  0518 

703  NSTFFM  CSCThCTAI  0519 

HSTPI-STN(THETA1  0520 

WF=MSTfF  C521 

NT  -  MS  ttt  0522  ' 

I  MI  N~  2  T523 

JVIN^NPALF  0524 

r G  TC  7?5  0525 

710  .1  =  JMTN  C526 

Ml  720  !=IMIS,NTOT,NF?  0527 

sump- (7atariii*oat«ri  jn/?.  352s 

SUMI- ICATAICI1 *OATAI (Jll/2,  0529 

r  IF  R  =  ( r  AT  AP  ( 1 1  -0  AT  AR  ( Jl  1  /  ?  .  G  5  30 

rir  I -(GAT  Aim -GAT  81  (Jll/2.  C5J1 

TEi>PRrw6*SUM*MI*r.IFR  0532 

TrMPI^6I*SG'M-NR«Cirp  0533 

P  AT  AF  1 1 1  s  SL'Kfi*  TEMFe  0534 

r.ATAIU»sOiri*TEFFI  0535 

PATAK (JI=SIMR-TEMFF  0536 

oatatiji--cifi«te(‘pi  0537 

720  J-J*6F2  0538 

IMIN-IMN*1  P539 

JM1N-JMIN-1  0540 

MTEMP-mR *MSTPI  0541 

MPrHR.HSTPF-Ml’MSIFI  Cc42 

WIMl*4STPF»MTtMF  "543 

72c  IF  ( IMF -JMTNt  710  ,7  70 ,740  0*44 

7  30  I«-(IFRwri731, 7  1,0,  731  0545 

711  GO  '15  I-IFTM,MOT.Nr?  0546 

7  ’ 5  PAIA  1(  II  -  -TATAK  IT  0r4  7 

740  Nr2 ' M "*NP2  0548 

MCT5NTCT*KTGT  0549 


10/. 


J-HI0T*1  0550 

TMAX=N1CT/2»1  0651 

7  46  IhJN-IHAX-KHALF  0652 

I-IHIK  "553 

CO  TO  755  0554 

760  r.ATAK  « J»  =  DA  T  ARTI1  GS55 

OATAlIJM-CATAItll  0656 

756  I=I»1  "'57 

j=j-i  c^se 

IF  U-IHAXI  756,760 ,760  O' 59 

760  OATAfclJ>  =  DATARIIf>.  IM-OATAI1IHTN1  0'60 

OATAHJ>  =  0.  (Vfci 

IFII-JI770,7£0 ,780  0'62 

765  OATARIJl  =  0ATAR'tI»  O' 63 

OATAI 7JI--OAIAITII  0664 

770  1=1-1  0665 

J  =  J-1  0566 

TF(I-IMN)  776,775,765  0567 

776  <'iATAR(.t)  =  OfTARtIKIM+ftATAl(TMlM  3568 

DATAI<J1=0.  0565 

INAX  =  IMK  0570 

00  TC  765  0571 

700  r.ATARm=fl»TAR«l»  ♦CATAI11  »  0672 

r.ATAI  (1»  -0  .  0573 

r,C  1C  r-U0  0571, 

r  0576 

C  COMPit |6  fi  KCAL  TRANSFORM  FOR  THE  IND,  3RD,  ETC.  CINF4SICK  9*  0576 

C  FOKJliCATE  FYHMETRIES.  0577 

C  5678 

80C  |F(NF-RFV-2l5ufl»9C6»805  O670 

806  DO  860  13=1 , NTOT , KF?  0680 

I?HAX=T3*NF?-NF1  2681 

nr.  SCO  12= 13,I2MAX,NF1  0682 

IhAX- I?»NP1-1  0663 

I  MI  N  =  12*  I1RNG  0681, 

JMBX:I3*IT»NP1-IKIK  0585 

IRIT2-T 3 1820  ,823,810  0586 

810  JMA  X  =  J8  AX ♦ KF  2  0587 

820  TFCTnlh-2M6u,850  ,*30  0588 

830  J.imav.NPO  '"589 

CO  860  MMMH4X  0590 

r.ATAR  (II  =  DATARTJt  0591 

DAI  AT  1 !> =-F ATATT Jl  05S? 

PAT  J=J-1  065? 

C'fl  .1  1MAX  0594 

no  860  I  -  I  6  IN , IHAX , NFC  ?6C6 

DAIAMI1 -DATBR'IJT  C696 

DATA  I  ( 1 1 - -r ATA U J>  0657 

F6C  J-J-KFO  3558 

C  0559 


.1 


F NT  f.f  tGO«=  ON  E#0  OIKENSION  0600 

C601 

KFC-NFl  0602 

0603 

KFFfV-K  C606 

PFTUSN  0^05 

ESC  0606 


DATE 

ILMEI 


